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FOREWORD 


In  January,  1983.  the  Air  Force  Geophysics  Laboratory  <AFGL> 
awarded  contract  number  F19628-83-C-0056  to  the  Space  Dynamics 
Laboratories  <SDL>  at  Utah  State  University  <U5U>  in  Logan,  Utah, 
under  which  SDL  assists  in  the  infrared  measurement  program. 

MAPSTAR  is  part  of  AFGL's  Laboratory  Simulation  program. 
The  overall  goal  of  the  program  is  to  create  a  controlled- 
environment  data  base  of  infrared  emissions  resulting  from 
molecular  interactions  in  order  to  more  accurately  characterize 
and  identify  the  sources  of  infrared  radiation  in  the  atmosphere. 
As  a  result,  our  ability  to  identify  and  track  man-caused 
disturbances  in  the  atmosphere  with  infrared  instruments  will  be 
greatly  enhanced.  /“ 

Laboratory  models  are  being  cross-checked  with  a'btuAl 
atmospheric  observations.  Based  on  these  findings,  nuclear  and 
atmospheric  predictive  models  will  be  refined  so  that  we  may 
assess  the  effects  of  mankind's  activities  on  atmospheric 
chemistry  and  more  accurately  identify  artificial  sources  of 
infrared  radiation,  such  as  rocket  plumes. 

AFGL  has  contracted  the  Laboratory  Simulation  Program  to 

SDL. 

MAPSTAR  is  part  of  a  ground-based  remote  sensing  program 
designed  to  define  infrared  structure  in  the  OH  region  of  the 
mesosphere.  Several  sensors  will  be  employed  in  the  MAPSTAR 
effort.  In  the  program's  first  stage,  a  high-resolution  ground- 
based  interferometer  equipped  with  a  high-straylight  rejection 
telescope  was  used  to  examine  the  spectral  characteristics  of 


iii 


airglow  from  sources  in  the  mesosphere.  Data  are  also  being 
collected  with  a  Doppler  radar  and  an  Isocon  optical  imager. 

The  interferometer-telescope  sensor  was  designed  to  gather 
quantitative  data  on  the  structure  of  airglow  originating  in  the 
mesosphere,  particularly  the  structure  of  alternating  bands  of 
airglow  attributed  to  "gravity  waves"  propagating  through  the 
atmosphere.  Data  gathered  in  June,  1983,  apparently  revealed 
wavelengths,  periods  and  phase  velocities  for  the  gravity  waves, 
as  well  as  rotational  temperatures  in  the  bright  and  dark  bands 
of  airglow,  providing  a  wealth  of  new  information  about  the 
chemical  and  physical  behavior  of  the  mesosphere. 


SUMMARY 


A  bright  OH  Meinel  airglow  structure  event  was  recorded  on 
June  15,  1983  from  Sacramento  Peak,  New  Mexico.  The  MAPSTAR 
interferometer  was  operated  in  conjunction  with  low-light-level 
infrared  imaging  isocon  camera  system  provided  by  the  University 
of  Southampton,  England.  The  camera  was  co-aligned  with  the 
telescope  to  provide  an  infrared  video  "eye"  for  the  interfero¬ 
meter. 

The  structures  were  measured  at  elevation  angles  near  the 
horizon.  Apparent  wavelengths,  periods  and  phase  velocities  of 
24+1  km,  14+1  minutes,  and  28+2  meters/second  respectively,  were 
calculated  for  the  recorded  structure.  The  interferometer  data 
show  intensity  modulations  of  20-40  percent  within  the  structure. 
A  mean  rotational  temperature  of  165  degrees  Kelvin  was  calcu¬ 
lated  using  the  interf erometer  spectral  data  and  temperature 
modulations  of  5  to  10  degrees  Kelvin  were  recorded  in  phase 
with  the  intensity  modulations. 

This  report  discusses  the  structure  and  performance  of  the 
MAPSTAR  interferometer  and  telescope,  as  well  as  the  results 
obtained  using  the  sensor.  Data  collected  shed  new  light  on  the 
phenomena  under  invest i gat  ion  and  will  provide  a  foundation  for 
studies  in  the  future. 
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CHAPTER  I 
INTRODUCTION 


The  mesosphere  is  the  interface  region  between  the 
earths ’  inner  and  outer  atmospheres .  Occur rences  at  the 
mesopause  include:  the  temperature  gradient  makes  a  sign 
change,  the  atmospheric  pressure,  density,  and  mean 
molecular  weight  all  have  an  inflection  point  in  their 
respective  curves  [Banks  and  Kockarts  19731.  The  atmosphere 
makes  a  transition  from  a  fluid  to  free  molecular  flow  in 
t'his  region  which  accounts  tor  these  changes.  The  unique 
properties  of  the  mesospheric  region  are  of  great  interest 
in  under stand i ng  the  middle  atmosphere  and  its  influence  on 
the  energy  budget  of  the  earth. 

The  80  to  100  km  region  (mesosphere)  is  a  difficult 
reg ion  to  study  because  the  altitude  is  too  low  for  direct 
satellite  observations  and  too  high  for  direct  balloon  or 
airplane  measurements.  Ground  based  studies  are  hampered  by 
the  intervening  atmosphere.  There  are  some  relatively 
transparent  atmospheric  "windows"  in  the  near  infrared, 
however.  Hydroxyl  radicals  <0H)  reside  in  this  mesospheric 
region  in  sufficient  concentrations  to  radiate  a  large 
amount  of  energy  at  the  red  and  near  infrared  wavelengths 
(Baker  19781. 
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The  excitation  o-f  OH  is  caused  by  various  solar  and 
chemical  processes.  The  excited  OH  radical  is  a  complex 
v 1 br at i ona 1 -rotat i ona 1  system  which  emits  radiation 
spectrally.  The  spectral  radiation  distribution  is  near 
Max  we  1  1  -  Bo  1  t  zmann  in  nature  [Baker  1978];  therefore,  the  OH 
spectral  radiation  can  be  measured,  the  Boltzmann 
distribution  determined,  and  a  rotational  temperature 
calculated  [Ware  1980].  This  temperature  can  then  be  used 
to  help  understand  the  chemistry  and  physics  o-f  the  entire 
region.  Banks  and  Kackarts  11973]  show  that  mesospheric 
temperatures  during  the  summer  months  at  a  mid-latitude  site 
can  be  expected  to  be  between  150  *K  and  190  *K. 

Recent  studies  have  shown  that  mesospheric  optical 
radiation  called  airglow  has,  at  times,  exhibited  some 
’wavel ike‘  structure  [Taylor  et  al.  1980].  These  waves  have 
been  studied  using  photographic  and  photometric  methods. 
The  object  of  this  study  was  to  develop  and  utilize  an 
instrument  to  provide  quantitative  data  o-f  OH  rotational 
temperature  and  intensity  variations  during  these  airglow 
structure  events  . 

Airglow  Structure  Measurement  Background 

The  atmospheric  airglow  layer  has  been  observed  and 
studied  -for  many  years  using  a  variety  o-f  methods.  The 
studies  conducted  during  the  decades  -from  1930  to  1970, 
however,  failed  to  recognize  the  nature  of  the  airglow 
structure  phenomenon.  Rayleigh  C19311  was  among  the  first 
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to  recognize  the  difference  between  airglow  and  aurora.  He 
referred  to  the  enhanced  airglow  as  “non-polar  aurorae.” 
Photographs  of  the  airglow  structure  were  presented  in  1952 
by  Hoffmeister  119523$  but  again,  neither  the  identity  nor 
the  source  of  the  airglow  structure  was  understood . 
Chamberlain  119613  briefly  outlines  the  historical  efforts 
in  airglow  studies  up  until  1961.  During  the  decade  of  the 
60s  the  techniques  of  photometry  were  perfected  and  most 
optical  atmospheric  study  efforts  were  centered  around  these 
methods. 

Kieffaber  119233  presented  photographic  evidence  in  1972 
of  apparent  airglow  “waves”  and  “cells’”  in  the  750-900  nm 
wavelength  region  using  infrared  film  and  a  35-mrn  camera. 
She  proposed  that  the  airglow  stripes  originated  from  a 
disturbance  in  the  OH  layer.  Again  in  late  1972,  Peterson 
and  Kieffaber  119733  recorded  more  occurrences  of  structure 
at  their  mid-latitude  site  near  Albuquerque,  New  Mexico. 
The  photographically  recorded  events  were  also  tracked  with 
infrared  photometers  <1.65  and  2.2  Pm)  and  shown  to  be 
moving  between  20  and  AO  meters/second. 

In  1975,  Crawford  et  al.  119753  flew  an  image 
intensified  isocan  television  system  on  board  NASA’s  CV990 
aircraft.  Peterson  and  Kieffaber  119753  observed  on  the 
same  flight  with  their  cameras  and  photometers.  Both  groups 
recorded  “cloud-1 ikt"  airglow  structures.  Again  using  35-mm 
cameras,  Moreels  and  Herse  119773  measured  extensive  OH 


airglow  structure  over  Europe. 


Their  findings  were  similar 
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to  those  of  Peterson  and  Kie-ffaber.  Uavcs  on  the  order  of 
40-km  spatial  wavelength  appeared  to  be  traveling  at 
horizontal  speeds  ■from  13  to  20  me  t  er  s  /  sec  ond  .  Peterson 
119791  was  able  to  record  numerous  occurrences  from  1975 
through  1978  with  some  of  the  events  being  enhanced  enough 
to  see  the  structure  with  the  naked  eye.  The  University  of 
Southampton  Atmospheric  Physics  group  [Taylor  et  al .  1980 
and  Taylor  1984]  recorded  many  structure  events  with  the 
image  isocon  television  cameras  from  1973  through  1983. 
Using  radiometric  techniques,  Huppi  and  Baker  11976]  also 
recorded  OH  intensity  variations.  Takeuchi  and  Misawa 
C1981]  record  some  short-period  waves  o-f  OH  intensity  and 
rotational  temperature  using  a  tilting  -filter  photometer. 
This  study  used  a  narrow  field  of  view,  fast  scan  rate 
instrument.  The  spectral  resolution  was  rather  coarse, 
however,  (unable  to  resolve  the  base  line  between  adjacent 
band  lines)  and  the  measurements  were  taken  without  the  aid 
of  any  photographic  or  video  equipment  making  it  difficult 
to  identify  what  was  being  observed.  Airglow  structure  was 
also  recorded  by  Peterson  and  Adams  119831  during  a  total 
lunar  eclipse  in  the  summer  of  1982.  In  this  lunar  eclipse 
study  extensive  use  was  made  of  photographic  equipment  as 
well  as  a  vertical  sounding  radar. 

During  the  decade  of  the  1970’s,  interferometric- 
spectroscopy  applied  to  middle  atmospheric  research  matured 
as  a  measurement  science.  Interf erogr ams  processed  using 
computer-based  fast  Fourier  transform  (FFT)  methods  yielded 
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high-resolution  OH  airglow  spectra  -from  which  intensity  and 
rotational  temperatures  could  be  extracted.  Baker  C19781 
presents  an  excellent  summary  of  the  studies  conducted  in 
this  area.  However,  because  the  instruments  used  had  wide 
fields  of  view  (therefore  integration  over  large  areas),  low 
throughput  (therefore  long  integration  times),  or  low 
spectral  resolution  the  small  spatial  airglow  variations 
were  not  able  to  be  spec t r oscop i ca 1 1 y  measured  at  high 
resolut ion . 

The  references  cited  and  many  others  have  recorded  OH 
airglow  structure  events.  However,  none  of  the  researchers 
were  able  to  provide  a  measure  of  high-resolution  spectral 
changes  and  therefore  calculate  the  differences  in  OH 
rotational  temperature  of  the  dark  and  bright  band  "waves." 
This  dissertation  gives  the  design,  development,  and 
operation  of  a  special  high-throughput,  narrow  field  of 
view,  fast  scan  interferometer-spectrometer  which  can 
spectrally,  spatially,  and  temporally  resolve  the  OH  airglow 
emission  structure. 


Airglow  Structure  and  the  Theory 
of  Atmospheric  Gravity  Waves 


In  a  landmark  paper,  Hines  C 19601  suggested 
certain  conditions  the  atmosphere  could  be  di 
“gravity  wave."  Later,  Hines  C19653  hypothesize 
passage  of  internal  gravity  waves  (IGW's) 
atmosphere  would  cause  some  reversible,  adiaba 
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(temperature  -fluctuations  associated  with  waves)  until  the 
dissipation  of  the  wave  became  excessive.  Along  with  his 
gravity  wave  hypothesis,  Hines  119653  gave  a  simple  model  to 
describe  his  suspected  IGW  temperature  fluctuations.  Using 
some  constants  from  the  h i gh -a 1 t i tude  vapor-trail 
measurements  of  Kochanski  119643  in  the  Hines  model,  a 
calculated  temperature  change  of  ±6  *K  could  be  expected  in 
conjunction  with  the  passage  of  an  IGW  through  the 
atmosphere.  Temperature  and  wind  measurements  by  Rai  and 
Fejer  119713  using  rocket -grenade  techniques  support  the  IGW 
hypothesis  put  forth  by  Hines. 

The  atmospheric  scientific  community  in  the  Soviet  Union 
has  done  extensive  work  in  the  area  of  an  adiabatic- 
oscillation  IGW  model  similar  to  the  work  of  Hines  C19603. 
Krassovsky  et  al.  119773  summarize  much  of  the  modeling  and 
measurement  efforts  of  this  group.  The  basic  instrument 
used  is  a  three-axis  diffraction  spectrograph,  supported  by 
various  photometers.  The  large  data  base  of  measurements 
presented,  shows  strong  correlation  between  the  magnitude  of 
the  temperature  modulation  and  the  period  of  the  IGW.  The 
observational  data  presented  exhibit,  for  IGW  periods  of 
about  25  minutes,  the  calculated  rotational  temperature 
changes  of  about  6  *K  with  the  passing  of  the  wave 
IKrassovsky  et  al.  19773.  Additionally,  the  Soviet  work 
indicates  that  a  significant  rotational  temperature 
difference  can  be  seen  between  high  level  and  low  level  OH 
vibrational  transitions.  The  observed  differences  range 
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■from  5  to  25  *K  with  the  high  vibrational  level  transitions 
appearing  hotter. 


Interferometers  and  Fourier 
Trans-farm  Spectroscopy 

The  Michelson  interferometer  was  invented  in  the  1880’s 
by  Albert  Abraham  Michelson  CShankland  19741.  In  the 
M iche 1  son -Mor 1 ey  experiment,  the  instrument  was  used  in  an 
attempt  to  measure  the  earth’s  movement  through  an  "ether. " 
Michelson  also  used  his  interferometer  to  determine  the 
exact  length  o-f  the  standard  meter  and  to  measure  the 
diameter  of  celestial  bodies.  He  also  discovered  the 
spectral  -fine  structure  o-f  hydrogen,  mercury,  and  thallium 
CShankland  19741.  This  pointed  out  the  potential  for  what 
would  be  later  be  called  "Fourier  spectroscopy." 

Opt i cal  1 y-sensi t i ve  detectors  were  subsequently  used  in 
conjunction  with  Michelson  interferometers,  in  which  one 
mirror  was  mechanically  displaced  at  a  constant  rate,  to 
produce  an  electrical  i nterf erogram.  The  interferogram  was 
then  inverted  using  Fourier  transform  techniques  to  yield 
direct  spectral  data.  In  1910,  Ruben  and  Woods  CConnes 
19631  obtained  the  first  far  infrared  spectrum  using  this 
method.  Fellgett  C19491  and  Jacquinot  C1954]  independently 
showed  the  inherent  advantage  that  the  interferometer  has 
over  grating  and  prism  spectrometers .  This  advantage 
results  from  measuring  all  spectral  components 
simultaneously.  This  improvement  was  referred  to  by 


Fellget  as  "multiplex  spectroscopy”.  Jacquinot  CConnes  19631 
also  showed  that  absence  of  slits  in  Michelson  spectrometry 
also  improved  system  throughput  when  compared  with  conventional 
grating  methods. 

The  advent  of  large,  fast  computers  and  the  development  of 
the  fast  Fourier  transform  CForman  19661 ,  which  together  could 
calculate  a  large  Fourier  transform  quickly,  spread  the  use  of 
Fourier  transforms  for  power  spectral  density  analysis  into 
many  fields.  Notable  contributions  were  made  in  the  field  in 
the  1950 ' s ,  1960 '  s ,  and  1970's  by  Kertz  £19591,  Connes  £19561, 
Gebbie  and  Vanasse  £19561 ,  Strong  and  Vanasse  £19591 ,  Forman 
£19661,  Haycock  and  Baker  £19751,  and  Steed  £19781.  Many  of 
the  improvements  in  the  field  were  reported  at  the  1970  Aspen 
Conference  on  Fourier  Spectroscopy  [Vanaaae  et  al .  19711. 

Figure  1-1  shows  the  typical  layout  of  the  Michelson 
interferometer.  The  system  is  comprised  of  a  beamsplitter  that 
divides  the  incoming  light  beam  into  two  equal  parts,  two 
mirrors  <M^  is  stationary  and  is  mobile),  and  a  condenser 
lens  to  focus  the  light  onto  a  detector. 

The  incoming  light  is  divided  into  two  portions  by  the 
beamsplitter  with  each  portion  directed  to  its  respective 
mirror.  The  energy  is  then  reflected  by  the  mirrors  and  is 
recombined  at  the  beamsplitter  and  passed  onto  the  detector 
system.  The  recombined  beam  is  modulated  by  differences  in 
the  path  lengths  between  the  beamsplitter  and  each  of  the 
mirrors.  If  the  mobile  mirror  M^  is  in  a  position  such 


MIRROR 


i 

INTERFEROGRAM 


Figure  1-1.  Layout  of  a  conventional  Michelson 
interferometer . 

that  the  path  length  AC  is  the  same  as  the  path  length  AB, 
then  the  recombined  signals  are  in  phase  and  thus  add 
construct ively .  The  same  constructive  addition  occurs  when 
the  path  length  difference  BC-AC  is  any  integral  number  of 
wavelengths  of  the  incoming  signal.  On  the  other  hand,  if 
the  path  length  difference  is  not  an  integral  multiple  of 
the  wavelength,  then  the  recombined  signal  will  have  varying 
amounts  of  destructive  interference  depending  upon  the  phase 
difference.  As  the  path  length  AC  is  changed  in  a  uniform 
manner,  by  moving  at  a  constant  rate,  the  electrical 
signal  from  the  detector  is  the  i nter f er ogram  of  the 
incoming  optical  signal.  The  Fourier  transformed 
i nterf erogram  yields  the  spectral  content  of  the  incoming 
1 i ght . 


The  simultaneous  measurement  o-f  high-resolution 

spectral,  temporal,  and  spatial  characteristics  o-f  the  OH 
airglow  structure  requires  an  instrument  with  both  a  narrow 
field  of  view  and  high  throughput.  The  standard  Michelson 
i nter f erometer ,  when  used  for  h i gh -reso 1 u t i on  measurements, 
has  a  narrow  field  of  view  but  its  low  throughput  would  make 
it  an  order  of  magnitude  less  sensitive. 

The  narrow  field  of  view  limitation  of  a  standard 
Michelson  interferometer  is  illustrated  in  Figure  1-2.  When 
incoming  energy  is  allowed  to  enter  the  interferometer  off- 
axis  (9#0*>  the  relationship  between  the  displacement  of 
mirror  M^  and  the  actual  path  difference  between  the  two 
mirrors  is  altered.  The  path  difference  or  retardation  is 
no  longer  2tf  as  it  is  when  light  is  coming  straight  into  the 
instrument,  but  now  is  2dcos0 ,  where  6  is  the  angle  of  the 
incoming  light  with  respect  to  the  entry  normal.  For  an 
instrument  with  a  given  resolving  power  the  maximum  usable 
field  of  view  for  a  standard  Michelson  i nter f erometer  is, 
according  to  Vanasse  C 19773 

n  =2w  t  R  ,  (1.1) 

max 

where 

A  =  maximum  field  of  view  in  steradians, 

max 

R  =  resolving  power  of  instrument. 

To  increase  the  throughput  (thus  achieving  a  faster  scan 
rate)  an  op t i ca 1 1 y -compensated  interferometer  was  chosen  for 
use  in  this  study.  There  have  been  many  proposed  methods 
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of  f i e 1 d -w i den i ng  or  optical-compensation  -for  Increasing  the 
throughput  o-f  an  interferometer.  These  techniques  are 
reviewed  by  Baker  CVanasse  19773.  The  method  used  -for  the 
instrument  in  this  study  was  -first  proposed  by  Connes 
C 19563,  and  uses  optical-compensation  wedges  or  prisms  in 
each  leg  o-f  the  interferometer  (specific  details  are 
discussed  in  Chapter  II).  Optical  compensation  increases 
the  throughput  by  increasing  the  usable  field  of  view  of  the 
instrument.  However,  the  measurement  of  OH  airglow 
structure  requires  a  small  field  of  view.  The  high 
throughput  of  the  compensated  i nterf erometer  was  matched 
(maintaining  temporal  resolution)  to  a  special  optical 
system  which  included  a  large-diameter  telescope  to  obtain 
the  desired  narrow  field  of  view  while  maintaining 
throughput . 

The  i nterf erometer  system  was  now  able  to  simultaneously 
resolve  the  spectral,  temporal,  and  spatial  characteristics 
of  the  OH  airglow  structure.  However,  because  the  measured 
radiation  was  in  the  infrared  it  was  necessary  to  locate  and 
measure  a  structure  occurrence  and  to  characterize  the  total 
structure  into  which  the  interferometer  was  looking.  The 
video  viewing  system  used  is  described  in  the  next  section. 


3 


Ieocon  Television  System 


Taylor  11983-843  participated  in  this  research  by 
providing  and  operating  a  low  light-level  infrared 
television  camera  used  in  conjunction  with  the 
interferometer-spectrometer.  This  Southampton  University 
infrared  TV  camera  allowed  any  airglow  structure  to  be 
quickly  and  efficiently  identified  and  permanently  recorded 
as  video  information.  Taylor's  system  employed  an  image 


intensified  isocon  television  system. 


This  i socon 


television  system  is  briefly  described  here  for 
comp leteness. 

The  isocon  tube  is  different  than  other  scanning  image 
tubes  in  that  it  uses  a  different  portion  of  the  scanning 
electron  beam  to  create  the  signal.  The  low-energy  electron 
beam  scans  a  h i gh-res i stance  target  as  in  other  tubes; 
however,  the  isocon  signal  comes  from  the  electrons  that  are 
scat tered  off  the  target  [Soule  19683  rather  than  the 
reflected  electrons  used  by  conventional  equipment 
(orthicons  use  the  reflected  beam).  The  scattered  signal, 
although  small  in  magnitude,  has  a  high  s i gna 1  -  to- no i se 
ratio,  and  is  particularly  well  suited  to  the  viewing  of  the 
low-contrast,  faint  OH  airglow. 

The  TV  system  used  to  image  the  airglow  structures  was 
specially  modified  to  enable  clear  images  of  the  OH  airglow 
patterns  to  be  obtained  in  less  than  a  second  [Taylor  19843. 
The  television  camera  used  was  an  English  Electric  Valve 
Miniature  Isocon,  Type  P1477  fitted  with  a  single-stage 


image  intensifier,  optically  coupled  to  a  55-mm  image  isocon 
tube.  The  camera  has  a  signal  -  to-noise  ratio  of 
approximately  40  dB  at  10  ft-candles  (starlight  conditions) 
and  a  dynamic  range  of  about  2000:1. 

To  further  enhance  the  capability  of  the  camera  to  image 
very  faint  airglow  structures,  it  was  arranged  for  the 
electronic  image  to  be  integrated  on  the  target  of  the  TO 
tube  for  a  period  of  up  to  a  second  (longer  integration 
times  allow  the  image  charge  on  the  target  to  migrate  and 
thus  smear  the  image)  before  being  scanned  and  recorded  onto 
video  tape  CTaylar  19843.  This  technique  is  particularly 
useful  as  it  improves  the  s i gna 1  - to- no i se  ratio  of  the 
airglow  signal  by  nearly  an  order  of  magnitude  (7  times  for 
a  1-second  integration  period)  with  no  significant  loss  of 
temporal  resolution. 

The  camera  has  an  extended  red  spectral  response  (S25) 
with  a  peak  sensitivity  at  about  500  nm  and  a  long 
wavelength  cut  off  at  900  nm.  Images  of  the  near-infrared 
OH  structure  were  obtained  by  placing  a  Schott  RG715  band 
stop  filter  in  front  of  the  camera  lens.  The  combined 
response  of  this  filter  and  the  TV  camera  gave  a  bandwidth 
(half  maximum)  of  715  to  850  nm  and  a  peak  sensitivity 
around  750  nm.  This  bandwidth  is  illustrated  in  Figure  1-3. 
The  location  of  all  the  OH  emission  bands  within  this 
spectral  range  are  indicated.  The  intensity  in  the  zenith 
of  the  OH  emission  within  this  bandpass  is  typically  5  to  15 
kR |  the  principal  emissions  are  the  OH  (9,4)  and  (5,1) 
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Figure  1-3.  Isocon  camera  system  spectral  range  including 
the  OH  emissions  within  that  range  tTaylor  1983-843. 

Meinel  bands  CTaylor  19843.  The  camera  was  fitted  with  a 
Nikon  85-mm,  f/1.4  lens  and  was  adjusted  to  give  an  almost 
square  field  of  view  of  15*  horizontal  by  13*  vertical.  The 
isocon  camera  used  in  this  study  is  shown  in  Figure  2-9, 


mounted  on  the  i nterf erometer  telescope. 
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The  Southampton  TV  cameras  have  been  used  tor  many  years 
to  "photograph"  the  near - i nf r ar ed  OH  airglow  structure.  An 
example  of  the  quality  of  the  video  data  gathered  is  shown 
in  Figure  1-4.  This  video  -frame  was  taken  in  August  of  1980 
while  observing  over  the  Swiss  Alps.  The  bright  and  dark 
bands  each  subtend  about  1*  of  arc.  These  data  were  the 
basis  upon  which  the  interferometer’s  field  of  view  was 
designed  in  order  to  resolve  the  spatial  nature  of  the  OH 
airglow  structure. 


Scope  and  Objectives 

The  specific  goals  and  objectives  are  outlined  as 
foil ows : 

1.  Design  and  develop  an  optical  instrumentation  system 
capable  of  quantifying  the  spatial,  spectral,  and 
temporal  characteristics  of  OH  near-infrared  night 
airglow  structure.  The  instrumental  field  of  view  must 
be  one  degree  or  less  to  resolve  the  structural 

characteristics  of  the  airglow.  The  system  must  have 
spectral  resolution  of  better  than  3  cm  1  in  order  to 
provide  spectra  from  which  OH  rotational  temperatures 
can  be  calculated  using  appropriate  algorithms  and 
digital  computer  programs.  The  NESR  must  be  sufficient 
to  resolve  the  OH  near - i nf r ar ed  airglow  with  scan  times 
of  less  than  1  minute  to  resolve  the  temporal 


fluctuations  of  the  airglow  structure. 


Figure  1-4.  Image  isocan  photo  o-f  OH  airg 
taken  by  Taylor  et  al.  [1980]  in  Switzerland  i 


Use  the  instrumentation  system  to  measure  the  spectral, 
spatial,  and  temporal  characteristics  of  OH  near- 
infrared  airglow  structure  -from  a  mid-latitude  observing 
site. 

Develop  and  apply  signal  processing  procedures  extending 
the  sampling  and  FFT  work  of  Ware  C19803  to  extract  both 
radiance  and  rotational  temperature  variations  of  OH 
airglow  structure. 

Derive  error  bounds  on  the  measurement  data,  based  upon 
system  specifications  such  as  field  of  view,  scan  speed, 
and  spectral  resolution,  as  well  as  instrument 
calibration,  and  signal  processing  techniques. 

Present  the  observational  results,  correlate  intensity 
variations,  temperature  fluctuations,  and  structure  with 
simultaneous  near-infrared  video  images,  and  compare  the 
findings  with  expected  OH  airglow  dynamics  studies  from 


other  i nvest i gator s . 


CHAPTER  II 


OPTICAL  INSTRUMENTATION  SYSTEM  DESIGN 

Design  Philosophy 

The  goal  of  this  study  was  to  develop  a  technique  -for 
simultaneously  measuring  the  spatial,  spectral,  and  temporal 
character  i  st  i  cs  o-f  OH  near- i  nf  rared  airglow  structures.  The 
basic  instrumental  approach  chosen  to  provide  the  spectral 
resolving  capability  is  a  Michelson  interferometer- 
spectrometer  which  is  optical ly-compensated  to  achieve  a 
very  high  throughput.  The  compensation  technique  used  makes 
it  possible  -for  obliquely  incident  optical  energy  up  to  5 
degrees  off  axis  to  contribute  to  the  detected  signal 
without  sacrificing  spectral  resolution.  This  resulting 
high  throughput ,  within  the  i nterferometer ,  is  matched  at 
the  input  to  a  large  diameter  collecting  telescope  yielding 
a  sub-degree  field  of  view  needed  for  spatial  resolving 
power.  The  entire  optical  system  is  diagrammatical ly  shown 
in  F i gure  2 - 1 . 

An  optically-compensated  interferometer  has  the  high 
throughput  needed  to  achieve  a  relatively  high  temporal 
resolving  power,  in  other  words,  a  scan  time  of  less  than  a 
minute.  Temporal  variations  are  therefore  identified  while 
maintaining  a  spectral  resolution  of  2  cm  *  in  the  near- 
infrared.  A  spectral  resolution  nearly  this  high  is 


Optical  layout  o-f  interferometer  telescope 
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desirable  in  order  to  unambiguously  compute  OH  rotational 
temperatures  from  the  measured  spectra.  A  free  spectral 
range  of  0.8  to  1.6  Mm  has  been  obtained  using  a 
cryogen ical ly-cooled  intrinsic  germanium  detector.  The 
optical  system  was  designed  to  maintain  the  high  throughput 
capability  (0.285  cm2  sr)  of  the  interferometer  while 
operating  at  a  narrow  field  of  view  (<1* )  in  order  to  be 
able  to  resolve  the  spatial  nature  of  the  airglow  structure. 

The  design  criteria  and  the  resulting  design  for  an 
1  optically-compensated  i nterf erometer  are  given  in  this 

i 

chapter.  Then  an  analysis  is  made  of  the  resulting 

i nterf erometer -spectrometer  system. 

i 

* 


;  High-Throughput  I nterf erometer  Design 

I 


The  optical  layout  of  a  conventional  Michelson 

i nterf erometer  was  discussed  in  Chapter  I.  Referring  to 

Figure  1-2,  the  retardation  or  path  difference  is  a  function 

of  the  entry  angle  of  the  incoming  energy.  The  relationship 

is,  Ac2</cos9,  where  o'  is  the  on-axis  (0=0*)  drive  distance 

and  0  is  the  angle  of  the  off-axis  ray.  The  maximum  field 

of  view  for  a  standard  Michelson  interferometer  is, 

ft  =2w / R  <Eq  1.1),  where  ft  is  the  maximum  field  of  view 

max  max 

of  the  instrument  and  /?  is  the  resolving  power.  When  viewing 
faint  airglow  events  this  limitation  on  throughput  is  a 
severe  one.  The  optical  system  is  compensated  to  increase 
system  throughput  for  temporal  resolution.  This  high- 
throughput  has  been  matched  to  a  collector  system  to 
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Figure  2-2.  Optical  compensation  method  conceived  by  Cannes 
£ 19563 . 


maintain  the  throughput  while  narrowing  the  -field  of  view 
(■for  spatial  resolution). 

The  method  used  in  this  study  to  improve  the  throughput 
of  the  interferometer  system  is  the  one  analyzed  by 
Bouchareine  and  Connes  119633  and  is  depicted  in  Figure  2-2. 
The  compensated-opt i cs  design  used  follows  that  developed  by 
Despain  et  al .  C19713,  and  the  design  limitations  previously 
derived  will  be  summarized  here  for  completeness.  Figure  2- 
3  shows  the  compensation  analysis  approach  of  Steed  £19783. 
With  the  optical  compensation  prisms  inserted,  the 
retardation  is 


23 


=  2d  cos  6  +  2 tn  cos  ♦  -  2d  cos  9  ,  (2.1) 


where 


4^=  retardation  with  compensation  prisms  inserted, 
t  =  thickness  of  optical  material, 
n  =  index  of  refraction  of  optical  material, 
d  -  mirror  drive  distance,  AC  -  AB, 

9  =  angle  of  off-axis  ray, 

♦  =  refraction  angle  of  off-axis  ray  through  optical 
mater i al . 


The  retardation  equation  shown  above  can  be  expanded  in 
a  Taylor  series  and  like  terms  collected  to  show  the  field 
of  view  dependency  more  directly;  namely, 


=  2 [d(/?  -  1)  +  tf]  ♦  [d 


n  —  1 


-  d]9‘ 


d  ^  1 1?  -  1  >  t  t  4 

-  +  -  +  -  -  -  9 

12  3n3  12j73  12 


(2.2) 


which  is  a  quartic  equation  in  9. 

Ideally,  the  retardation  A^_  should  be  independent  of  the 
entry  angle  9.  Equation  2.2  shows  that  this  is  not  possible 


by  varying  only  the  drive  distance  d  and  the  index  of 
refraction  n.  However,  a  significant  improvement  can  be 
realized  by  designing  such  that, 


d  =  t  [ 


77  —  1  • 


(2.3) 


which  eliminates  the  9  term  in  the  retardation  Equation  2.2 


by  forcing  the  9  coefficient  to  zero. 


This  leaves  only 


V1. 


Figure  2-3.  Optical  compensation  with  an  optical  section 
inserted  into  one  leg  of  a  Michel  son  interferometer  CSteed 
19781  . 

I 

4 

the  G  term  and  since  0  is  small  this  term  is  relatively 
i  ns i gn i f i can t . 

Designing  the  compensation  wedges  to  satisfy  Eq  2.3  and 

then  solving  Equation  2.2  for  an  on-axis  ray  (6=0*),  the 

resulting  retardation  A  is 
l  co 

(2.4) 

front  and  using  the 
retardation  for  the 


Assuming  a  monochromatic  wave 
results  of  Eq  2.4,  the  general 


compensated  case  is 
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where 

A  =  compensated  retardation, 

A  =  retardation  for  an  on-axis  ray. 

0  =  entry  angle  tor  incoming  ray  in  radians 

r>  =  index  of  retraction  ot  compensation  prisms. 

The  analysis  tor  the  Connes  [19563  method  shown  in 
Figure  2-2,  in  which  only  one  optical  component  is  driven, 
is  the  same  as  tor  the  system  shown  in  Figure  2.3.  The 
thickness  ot  the  optical  material  must,  however,  increase 
with  the  drive  distance  d  in  order  to  maintain  optical 
compensation  (see  Equation  2.3).  The  increase  in  optical- 
material  thickness  needed  to  maintain  compensation  is 
obtained  by  driving  one  ot  the  optical  wedge/mirror 
assemblies  in  synchronism  with  the  drive  motor.  Reterring 
to  Figure  2-2,  it  can  be  seen  that  the  reflective  elements 
or  mirrors  in  this  method  are  created  by  depositing  the 
mirrored  surface  to  the  back  side  of  each  wedge. 

Steed  [19783  derived  the  limits  on  field  of  view  as  a 
function  of  resolving  power  based  upon  aberration  limits. 
The  results  are  categorized  into  several  groups.  The  first 
is  chromatic  aberration,  that  is,  differences  in  the 
compensation  because  the  index  of  refraction  of  the  optical 
material  is  a  function  of  wavelength.  The  chromatic  limit 
i  s 


fl  „  n  ( n  -  1  )  x 


A  r  “m  "  1 

c  I  5  n  a 

L  co 


(2.6) 
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where 

fl 

( 

ft 


M 


*  field  of  view  at  wavelength  X  in  steradians, 

=  field  of  view  for  conventional  Michelson 
i nter  f erometer , 

n  =  index  of  refraction  at  wavelength  X, 

Sn  =  rt-Tt  where  n  =  index  of  refraction  at  compensated 
c  c 

wave  1 ength , 

#  =  desired  resolving  power, 

A  =  retardation  for  on-axis  ray, 

co 

X  =  wavelength  expressed  in  same  units  as  A 

The  second  limiting  factor  on  field  of  view  is  that  of 

spherical  aberration.  This  limit  is  also  derived  by  Steed 
[19781.  The  maximum  field  of  view  assuming  only  spherical 
aberrations  is 


ft  =  11  n  J  (2/?) 
s  M 


(2.7) 


The 

third  type 

of  aberration  is 

ast i gmat i sm. 

In 

the 

Connes 

method  a 

wedge  is  placed 

in  each  leg 

of 

the 

interferometer;  however,  only  one  of  the  wedges  is  driven. 
This  simplifies  the  mechanical  design  but  because  a  wedge  is 
in  each  optical  path,  astigmatism  aberrations  occur.  The 
severity  of  the  aberration  increases  as  the  angle  (cc )  of  the 
compensation  wedges  increases.  Bouchareine  and  Connes 
C 1963 3  quantify  this  aberration  as 

(1  *  2ft  /tan2a  ,  (2.8) 
A  M 
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where 

8  =  prism  angle  of  the  compensation  wedge. 

The  aberrations  considered  were  each  derived  assuming 
the  net  effect  o-f  the  distortion  was  a  shift  of  one  fringe 
in  the  interference  pattern.  These  degradations  prove  to  be 
the  limiting  factors  upon  the  maximum  usable  field  of  view. 
Figure  2-4  is  a  plot  of  Equations  2.6,  2.7,  and  2.8  and 
shows  the  relationship  between  these  various  limits.  At 
least  an  order  of  magnitude  improvement  in  throughput  is 
obtainable  using  the  wedge  prism  compensation  technique. 

The  next  step  in  the  system  design  is  to  ascertain  the 
optical  retardation  for  the  Corines  119561  method.  Referring 
to  Figure  2-5,  Steed  C1978D  showed  the  drive  plane  is 
parallel  with  the  apparent  image  plane  of  the  optical 
wedges.  The  difference  in  optical  path  length  or 
retardation  between  image  points  A1  and  A2  separated  by  a 
drive  distance  X  is 

6  =  2X  sinip-y)  ,  (2.9) 

where 

P  =  sin  *(;?  sina), 

Y  =  8  -  tan  A[<- - ->  tana], 

L  n  i 

a  =  wedge  angle, 

X  =  drive  distance. 


FIELD  OF  VIEW  (Steradians) 


RESOLVING  POWER 

Figure  2-4.  Plot  o-f  aberration  limits  -for  Connes  C1956 
method  o-f  optical  compensation  compared  with  conventiona 
Michelson  interferometer  [Steed  1978]. 


Drive  System 


Examination  of  Equation  2.9  shows  that  the  choice  of  a 
large  wedge  angle  reduces  the  actual  drive  distance,  thereby 
easing  the  mechanical  drive  requirements.  However,  the 
choice  must  be  made  in  conjunction  with  the  limits  on 
distortion  presented  in  Figure  2-4.  The  wedges  and  cube 
beamsplitter  used  in  this  instrument  are  made  of  a  high 
quality  quartz,  namely,  Infrasil  I,  manufactured  by  Amersil 
Corporation.  The  measured  index  of  refraction  for  this 
material  is  1.45.  A  compensation  prism  wedge  angle  of  8* 
was  chosen.  Solving  Equation  2.9  far  this  wedge  angle,  then 
p  =  11.64*,  Y  =  5.50°,  and  A  =  0.21X.  The  optical 
retardation  for  this  compensation  technique  is  only  about 
20%  of  the  drive  distance,  where  as  in  a  conventional 
Michelson  i nter f erometer  the  retardation  is  equal  to  2d 
where  d  is  the  drive  distance.  Thus,  for  a  given  spectral 
resolution  the  compensated  interferometer  requires  a  drive 
distance  of  2/. 21=9. 5  times  longer  than  the  conventional 
Michelson  approach. 

An  optically-compensated  i nterf erometer  with  a  spectral 
resolution  of  2  cm  1  requires  a  slide  movement  of  3.1  cm. 
This  comparatively  large  drive  distance  is  accomplished 
using  a  gas- lubr icated  platform,  providing  near-zero 
friction,  developed  by  Haycock  C19751.  The  platform 
translates  the  optical  wedge/mirror  assembly  and  as  such 
must  maintain  mechanical  tolerances  as  close  as  possible  to 
the  wavelength  dimensions  of  the  infrared  light  being 


measured.  The  gas  bearing  surfaces  are  lapped  to  a 

tolerance  of  about  1  wavelength  of  5461  A  light  and  operate 

-4 

with  a  clearance  of  2.5x10  cm.  At  an  operating  pressure  of 
6-12  psi  the  bearing  will  maintain  an  optical  alignment  of  1 
arc  second  while  translating  5  cm  [Haycock  19753.  The 
optical  components  are  then  mounted  on  the  bearing-supported 
optical  platform.  Adjustment  for  parallelism  is 

accomplished  on  an  optical  bench  using  a  laboratory  He-Ne 
laser  as  a  light  source.  The  material  used  in  the  platform 
and  bearing  assembly  is  a  specially  formulated  Invar  alloy, 
chosen  for  its  excellent  temperature  and  stability 
charac ter i st i cs  and  to  match  the  thermal  expansion 
character i st i cs  of  the  optical  material.  This  sophisticated 
platform  provides  the  relatively  long  drive  distances 
required  (up  to  5  cm)  as  well  as  the  mechanical  accuracy 
needed  to  translate  the  optical  components  properly. 

The  platform  is  driven  by  a  “voice  coil"  motor  and 
therefore  makes  no  physical  contact  with  the  rest  of  the 
instrument.  The  motor  is  driven  by  a  standard  servo- 
amplifier  with  feedback  from  positional  and  velocity  sensors 
located  within  the  platform  slide.  The  main  signal  driving 
the  slide  is  generated  by  a  digitally-controlled  ramp  which 
can  be  adjusted  for  slide  velocity  and  drive  distance.  The 
slide  controller  can  produce  variable  scan  times  of  5 
seconds  up  to  several  minutes,  and  a  drive  distance  of  up  to 
5  cm.  The  completed  i nterf erometer  uses  a  10-cm  cube 
beamsplitter  and  11.4-crn  diameter  end  wedges  with  a  prism 
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angle  of  8*.  The  interferometer  with  bearing  system  is 
shown  in  Figure  2-6. 

For  simplicity,  a  secondary  HeNe  1 aser -exc l ted 
interferometer  is  used,  with  its  independent  moving  mirror 
on  the  same  optical  plat-form  as  the  primary  optical  signal 
channel,  to  monitor  the  slide  position.  The  laser  signal  is 
counted  down  by  6  and  used  to  digitize  the  main  channel 
i nterf erogram.  This  method  provides  sufficient  sample 

points  for  a  16k  fast  Fourier  transform. 

The  other  spec i f i cat  ions  for  accuracy  of  the  optical  and 
mechanical  components  are  described  in  detail  by  Steed 
C19783.  The  entire  system  was  originally  designed  to  be 
operated  at  liquid  nitrogen  temperature  (77  »K)  to  reduce 

background  radiation  at  longer  wavelengths  ( >2  Mm).  The 
spectral  range  used  in  this  study  (0.8  to  1.6  Mm)  is  not 
background  limited  and  therefore  does  not  require  the  added 
complexity  of  optical  train  cooling. 

Detector  System 

The  optically-compensated  i nter f erometer  has  a  very 

2 

large  throughput  (Ad  «  l  cm  sr).  The  desired  goal  of 

optical  compensation  was  to  improve  the  sensitivity  of  the 

i nter f erometer  system!  however,  to  take  advantage  of  the 

large  throughput  a  large  diameter  detector  is  required.  For 

2 

example,  assuming  a  throughput  of  1  cm  sr  and  collector 
optics  with  a  very  fast  f-number  of  0.5,  the  detector 
diameter  would  have  to  be  0.650  cm.  Large  infrared 
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detectors  with  the  necessary  detectivity  are  very  difficult 
to  manufacture;  therefore,  some  compromise  was  necessitated. 
The  required  detector  diameter  and  the  speed  of  the 
collecting  optics  shown  in  the  example  are  both  somewhat 
impract ical . 

A  detector  system  was  chosen  which  was  readily  available 
and  had  proven  effective  in  the  past.  The  detector  selected 
was  an  RCA  Ltd.,  solid  state  germanium  device  with  a  self- 
contained  preamplifier  and  load  resistor  in  a  liquid- 
nitrogen  dewar.  Figure  2-7  is  a  picture  of  the  detector 
dewar .  The  detector  is  5  mm  in  diameter  and  has  a  noise 
equivalent  power  (NEP)  of  1.12  x  10  W/fHz  at  1.27  pm  and 
covers  a  spectral  range  of  0.8  to  1.6  pm. 

It  was  decided  that  a  scan  rate  of  30  seconds  would  be 
sufficient  to  resolve  the  temporal  variations  of  first 
interest  in  the  OH  airglow  structure.  This  decision  was 
based  upon  video  data  gathered  by  Taylor  et  al .  C 19801.  A 
30-secand  scan,  a  minimum  wavelength  of  0.8  Pm,  and  a  drive 

i 

distance  of  3.1  cm  (spectral  resolution  of  2  cm  *)  set  the 

i 

|  detector  electrical  bandwidth  at  150  Hz.  The  RCA  detector 

l 

has  a  bandwidth  of  600  Hz.  Appendix  A  contains  the  detailed 

| 

I 

’  specifications  of  the  detector  system. 

i 

! 

i  Telescope  Design 

I 

1 

i 

The  spatial  char ac ter i st i cs  of  the  airglow  structure, 
shown  by  Taylor  11983-841,  dictated  that  the  instrument 
field  of  view  be  less  than  1*  full  angle.  The  goal  of  the 
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optical  design  was  to  maintain  the  inherent  high  throughput 
of  the  compensated  i nterf erometer  throughout  the  entire 
system.  All  portions  of  the  optical  path  were  examined  in 
order  to  ascertain  where  the  system  was  throughput  -  1 i m i ted . 
The  choice  of  a  detector,  and  practical  limitations  on  the 
f-number  of  available  condenser  lenses  proved  to  be  the 
limiting  factors.  An  optical  system  was  designed, 
incorporating  the  interferometer,  using  a  large  diameter 
collector,  which  maintained  the  system  throughput  while 
narrowing  the  field  of  view  to  under  1*.  The  instrument 
system  was  now  capable  of  measuring  the  temporal,  spectral, 
and  spatial  variations  of  the  airglow  layer. 

The  design  of  the  optical  system  began  with  the  detector 

and  choice  of  its  condenser  lens  and  proceeded  back  towards 

the  telescope  collector.  The  entire  optical  system  Is 

described  d i agr ammat i ca 1 1 y  in  Figure  2-1.  The  diameter  of 

the  condenser  lens  must,  however,  allow  the  converging  beam 

to  pass  through  the  beamsplitter  and  wedges  without 

vignetting.  A  commercially  available  condenser  lens  with  a 

diameter  of  68  mm  and  an  effective  focal  length  of  50  mm  was 

chosen.  The  condenser  lens  and  detector  diameter  of  5  mm 

2 

set  the  system  throughput  or  An  at  0.28  cm  sr.  Referring  to 

2 

Figure  2-1,  a  throughput  of  0.28  cm  sr  translates  to  a  beam 
diameter  entering  the  interferometer  of  92  mm  and  converging 
in  a  6*  full  field  of  view.  The  physical  dimensions  of  the 
detector  optics  outlined  here  fit  well  within  the  size 


constraints  of  the  basic  interferometer.  The  field  of  view 
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is  also  within  the  aberration  limits,  shown  in  Figure  2-4, 
•for  an  optically-compensated  interferometer. 

The  calculated  throughput  can  now  be  used  to  design  the 

telescope  needed  to  narrow  the  field  of  view  to  1*  or  less. 

2 

A  throughput,  dictated  by  the  detector  system,  of  0.28  cm  sr 
and  a  1*  field  of  view  design  goal,  set  the  collector 
diameter  at  D  =  Att/2v  (l-cos6>=44  cm.  A  50.8-cm  (20-inch) 

diameter  system  was  chosen  because  of  availability  and  to 
allow  for  some  error  in  aligning  the  telescope  optics. 
Taking  advantage  of  the  larger  collector,  the  field  of  view 
was  narrowed  to  0.8*  requiring  a  48-cm  diameter  collector, 
within  the  50.8-cm  mirror  size  and  still  allowing  for  some 
error  in  optical  system  alignment.  For  portability,  an  f/2 
primary  was  specified  to  minimize  the  telescope  length.  An 
overall  system  f  number  of  8  was  selected  to  transfer  the 
telescope  image  to  the  proper  position. 

In  this  application,  the  purpose  of  the  telescope  is  to 
gather  energy  and  transfer  it  to  the  detector,  rather  than 
to  transfer  a  spatial  image.  Therefore,  the  quality  of  the 
optical  image  within  the  field  of  view  is  not  of  as  great  a 
concern  as  it  would  be  in  an  imaging  system.  Many  types  of 
folded  telescope  systems  were  considered.  In  order  to 

simplify  the  optics  Driscoll  and  Vaughan  [19731  suggests  a 
Dal 1-Kirkham  type  because  of  the  spherical  secondary,  if  the 
resulting  image  distortion  is  acceptable.  The  distortion 
for  a  Da  1 1 -K 1 r Wham  system  of  this  size  was  calculated  using 


the  formula  in  Driscoll  and  Vaughan  [19731,  and  found  to  be 
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less  than  0.4%  of  the  area  (coma  and  astigmatism  distortion 
were  calculated  in  terms  of  primary  mirror  area  that  will  be 
degraded)  o  f  the  primary  mirror.  Therefore,  a  Dali  -Kirkham 
type  o f  telescope,  being  more  than  adequate,  was  chosen  -for 
its  simplicity  and  relatively  low  cost.  This  type  of 
telescope  has  a  spherical  secondary  mirror  and  an  elliptical 
primary  mirror.  The  telescope-equipped  interferometer 
system  with  the  isocon  camera  mounted  on  the  telescope  is 
pictured  in  Figure  2-9. 

The  last  step  in  the  optical  design  was  to  provide  an 
optical  interlace  between  the  telescope  and  the  detector 
sub-systems.  Two  primary  cons i der at i ons  were  given  emphasis 
in  this  interlace  design:  (1)  imaging  the  detector  on  the 
primary  mirror,  and  (2)  imaging  the  condenser  lens  on  the 
telescope  Focal  plane.  The  detector  may  have  an  uneven 
response  across  its  area,  therefore  imaging  the  detector  on 
the  primary  mirror  minimizes  the  effects  of  off-axis  rays  by 
illuminating  the  entire  detector  by  light  entering  within 
the  field  of  view  of  the  instrument.  An  image  of  the 
detector  condenser  lens  at  the  telescope  focal  plane  allows 
for  independent  control  of  the  system  field  of  view. 

An  adjustable  iris  was  then  placed  at  the  focal  plane 
for  adjusting  the  field  of  view  while  still  allowing  the 
entire  detector  to  be  illuminated  by  energy  within  the  field 
of  view.  Careful  choice  of  physical  dimensions  and  optics 
allow  the  transfer  of  each  image  independent  of  the  other. 
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The  goal  of  placing  a  detector  image  on  the  primary 
mirror  was  accomplished  by  the  use  of  two  lenses.  Referring 
to  Figure  2-1,  a  transfer  lens,  which  has  a  diameter  o-f  100 
mm  and  a  -focal  length  o-f  250  mm,  is  placed  near  the 
interferometer  beamsplitter.  This  lens  produces  an  image  o-f 
the  detector  at  a  focal  point  near  the  physical  entrance  to 
the  i nterf erometer  housing.  The  magnification  ratio  of  this 
transfer  lens  focal  length  to  the  focal  length  (50  mm)  of 
the  condenser  lens  determines  the  detector  image  size  at 
this  point;  therefore,  the  detector  image  is  25  mm  in 
diameter.  A  field  stop  is  placed  at  this  point  to  limit  the 
detector  image  size  throughout  the  rest  of  the  system. 
Another  transfer  lens,  which  has  a  75-mm  diameter  and  a 
focal  length  of  200  mm,  is  placed  at  the  focal  plane  of  the 
telescope.  This  lens  transfers  the  detector  image  at  the 
field  stop  onto  the  primary  mirror.  As  can  be  seen  from 
Figure  2-1,  two  angled  mirrors  were  required  to  meet  the 
physical  constraints  of  the  design. 

The  second  design  abjective  of  the  optical  transfer 
system  was  to  place  an  image  of  the  detector  condenser  lens 
on  the  telescope  focal  plane.  A  collimated  image  of  the 
condenser  lens  is  generated  by  the  lOO-mm  diameter  transfer 
lens  that  is  placed  near  the  i nter f erometer  beamsplitter. 
The  addition  of  a  75-mm  diameter,  200-mm  focal  length  lens 
placed  at  the  field  stop  focuses  this  image  at  the  telescope 
focal  plane.  The  placement  of  this  lens  at  a  detector  image 
point  does  not  affect  that  image.  An  adjustable  iris  was 
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then  placed  at  the  telescope  -focal  plane  -for  control  of  the 
field  of  view.  The  field  stop  and  the  iris  provide  the 
limiting  apertures  for  the  entire  system.  The  optical 
system  described  here  was  verified  by  Harris  [19843  using  a 
computer-aided  optical  ray-tracing  program. 

The  optical  lenses  used  in  the  design  are  made  of 
commercial-grade  optical  glass  and  have  no  coatings.  The 
beamsplitter  and  wedges  are  constructed  of  quartz  and  are 
not  coated.  The  primary  and  secondary  telescope  mirrors  are 
aluminized  reflection  surfaces  with  a  S i 0^  coating. 

Instrument  Housing 

The  interferometer-spectrometer  with  its  associated 
optical  and  telescope  systems  was  placed  in  a  30-inch 
diameter,  45-inch  tall,  round  container.  The  interferometer 
placed  in  the  instrument  housing  is  shown  in  Figure  2-8. 
The  telescope  attached  to  the  container  is  shown  in  Figure 
2-9.  The  large  container  size  was  chosen  to  support  the 
telescope  in  a  stable  manner.  The  base  of  the  package  was 
mounted  upon  a  rotating  stand  which  provides  the  ability  to 
move  the  telescope  in  azimuth.  The  rotating  base  has 
attached  three  large  pneumatic  tires  for  system  mobility  and 
three  leveling  screws  to  stabilize  the  i nterf erometer  for 
operation.  The  telescope  is  mounted  to  the  side  of  the 
container  with  a  lO-inch  diameter  ball-bearing  allowing  the 
telescope  to  be  rotated  in  elevation.  The  two  rotating 
joints  provide  the  telescopy  with  complete  pointing  freedom. 
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The  specifications  of  the 
are  given  in  Table  2-1. 


instrument  with  the  telescope 


TABLE  2-1.  Summary  of  i nterf erometer-spectrometer 

spec  if i cat  ions. 


Throughput  ........ 

Scan  period  (minimum)  .  . 

Collector  diameter  .  .  .  . 

Field  of  view  (full  angle) 

Spectral  range  . 

Spectral  resolution  .  .  . 

Detector  type  . 

Detector  NEP . . 

Dynamic  range  . 

System  sensitivity  .  .  .  . 


^  2 

.  0 . 28  cm  sr 

.  1  scan/  30  seconds 

.  50 . 8  cm 

.  0 . 8  degrees 

.  0.8  to  1.6  Pm 

. 2  cm 

.  Intrinsic  germanium 

1.1  x  10  14  watts/fhz  at  1.5  pm 

. 80  dB 

.  16  R/cm  1  at  1.5  pm 


Figure  2  -  3  . 


Ir.tcrferameter  placed 


inside  housing. 


Figure  2-9.  Interferometer  system  equipped  with  20-ir.ch 
diameter  telescope.  The  isocon  camera  is  also  mounted  on 
the  tel  esc ope . 
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CHAPTER  III 


MEASUREMENT  THEORY 


The  I n ter ferogram 


The  electronic  signal  -from  the  i  nterf  erometer  detector 
is  a  law-frequency  time-continuous  signal  called  an 
i  nter-f  erograrn .  The  i  nterf  erogram  is  a  scaled  analog  o-f  the 
incoming  light  -frequency.  Using  the  approach  o-f  Loewenstein 
C19713,  and  referring  to  the  layout  of  a  Michelson 
i nterf erometer  in  Figure  1-1,  the  i nterf erogram  signal  is  of 
the  form 


f  =  t /<S  ,  (3.1) 

s 

where 

ff  «  input  light  wavenumber  in  cm 

v  =  rate  of  change  of  effective  path  difference  in 
cm/second , 

f  ■  scaled  frequency  output  of  interferometer  (the 

i nterf erogram)  in  Hz. 

The  detector  used  in  this  study  is  sensitive  in  the  region 
from  12500  to  6250  cm  1  <X  =  0.8  to  1.6  Mm).  The  scan  rate 

for  the  interferometer  was  set  at  30  seconds!  therefore,  the 
input  light  frequency  was  scaled  to  audio  frequencies  of 


less  than  200  Hz. 
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The  light  power  reaching  the  detector  in  an  ideal 
interferometer  is  ILoewenstein  19711 


P  ,  .  =  2A  tl+cos  2v ff x ) 

de  t 


(3.2) 


where 

A  =  amplitude  of  on-axis  monochromatic  point  source, 

x  =  effective  path  difference  in  cm. 

In  the  Loewenstein  model,  the  single  line  o-f  Equation  3.2  is 

replaced  by  an  input  power  density  spectrum  o-f  the  -form 
2 

A  =B  (ff ) .  Making  this  substitution  in  Equation  3.2, 

neglecting  the  constant  <dc)  part,  and  integrating  over  ff 
produc  es 


I  (x) 


B  ( ff  )  cos  (  2w  ff  x  )  tfff 
« 


(3.3) 


where  I ( x )  is  defined  as  the  i n t erf er ogr am .  The  desired 
information  B < ff )  (the  input  spectrum)  is  obtained  by  taking 
the  inverse  Fourier  transform  of  I(x).  The  signal  produced 
by  the  interferometer  is  referred  to  as  a  ’double-sided*  or 
’symmetric"  i nterf erogram.  An  example  of  a  double-sided 
i nterf erogram  is  shown  in  Figure  3-15  note  that  the  same 
information  is  available  on  both  sides  of  the  large  center 
which  occurs  at  zero  path  difference.  The  interferometer 
slide  is  moved  an  equal  distance  on  each  side  of  the  optical 
zero  path  difference  to  create  a  symmetrical  signal. 

Since  the  i n ter f erogram  I(x)  in  Equation  3.3  is 
symmetric,  the  spectrum  B(ff)  may  be  obtained  using  a  simple 
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Figure  3-1  Typical  "double  sided"  i  n  ter-f  erogram . 

cosine  transform  of  the  form 

B(ff)  =  |  I  <  x )  cos(2ifffx)  d  x  .  (3.4) 

*  • 

The  power  spectral  density  is  computed  with  a  digital 
computer  using  a  Fast  Fourier  Transform  (FFT).  The  FFT  used 
to  process  the  data  for  this  study  was  developed  by  Ware 
C1980I.  The  digitized  interf erogram  was  produced  using  a  12 
bit  analog  to  digital  converter  on  the  detector  amplifier 
output  and  enough  samples  were  taken  to  perform  a  16,384- 
point  transform.  The  implementation  of  the  FFT  is  beyond 
the  scope  of  this  paper  but  that  used  in  Fourier  transform 
spectroscopy  was  adapted  by  Forman  C 19661  from  the  radar 
signal  processing  work  of  Thomas  Stockham,  then  at  MIT 
Lincoln  Laboratory.  The  FFT  is  explained  extensively  by 
Brigham  C19743. 

The  computation  of  the  FFT  requires  that  the 
i nterferogr am  be  sampled  at  uniform  increments  of  path 
difference.  Ideally,  the  laser  reference  channel  should  be 
a  short  wavelength  laser  beam  which  passes  through  the 
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interferometer  optical  train  as  does  the  signal.  Otherwise, 
the  reference  beam  may  experience  different  motion  stability 
than  the  signal  beam,  due  to  nonidentical  geometry. 
However,  tor  simplicity,  a  separate  laser  (X  =  6328  A) 
i n ter t er ome ter  channel  was  used  and  provides  the  uniform 
sampling  tunction  signal.  The  laser  channel  signal  is 
divided  to  obtain  a  sample  rate  of  630  Hz.  The 
interf erogram  is  digitized  at  a  free  running  rate  of  30  kHz. 
The  values  taken  between  laser  channel  zero  crossings  are 
then  averaged  and  stored  as  the  data  point  for  the  slide 
position  cor respond i ng  halfway  between  adjacent  laser 
crossings.  This  over-sampling  method  will,  according  to 
Ware  119801,  minimize  the  system  noise  gain  as  well  as 
minimize  the  effects  of  slide  velocity  variations  in  the 
sampled  data. 


Calibration  Source 

The  determination  of  relative  OH  spectral  line  strengths 
is  sufficient  to  compute  rotational  temperatures.  Since 
this  is  a  major  interest  in  this  study,  a  calibration 
technique  to  produce  a  relative  instrument  response  was 
developed.  The  absolute  calibration  of  an  Interf erometer- 
spectrometer  is  an  intricate  process!  details  of  the 
approach  are  outlined  by  Wyatt  [ 19781.  Ware  119801 
performed  an  i n t er f er ome t er  wavenumber  response  calibration 
as  a  function  of  optical  alignment.  He  found  that  the 
relative  response  must  be  reestablished  each  time  the 
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instrument  is  realigned  at  the  time  of  data  collection 
because  the  instrument  Mould  only  remain  in  acceptable 
alignment  -for  about  2  hours.. 

A  relative  instrument  response  Mas  obtained  by  causing 
the  interferometer  to  view  a  calibration  source  composed  o-f 
a  tungsten  bulb  illuminating  a  pair  of  ground  glass 
diffusion  screens  (see  Figure  3-2).  The  intensity  of  the 
illumination  was  controlled  by  placing  a  plate  with  a  small 
3-mm  aperture  between  the  bulb  and  the  diffusion  screens. 
The  second  screen  was  then  imaged  by  a  projector  lens  onto  a 
24-inch  diameter  ground  glass  viewing  screen  placed  ten  feet 
away.  The  net  brightness  of  the  source  was  adjusted  so  the 
calibration  source  intensity  as  viewed  by  the  interferometer 
was  approximately  that  of  the  night  sky.  The  calibration 
source  assembly  was  mounted  inside  of  a  24-inch  diameter, 
lO-foot  long  tube.  The  viewing  screen  was  sufficiently 
large  to  completely  fill  the  interferometer  field  of  view, 
thus  providing  a  diffuse  source. 

The  long  length  (3.0  meters)  of  the  calibration  source 
was  chosen  to  ensure  that  the  illumination  on  the  viewing 
screen  was  uniform  across  the  entire  surface  viewed  by  the 
i nterf erometer .  The  quality  of  the  calibration  source  was 
calculated  using  the  geometry  described  above  and  assuming 
the  tungsten  bulb  was  a  blackbody  source  with  correction  for 


the  emissivity  of  tungsten. 
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VIEWING  SCREEN 


Figure  3-2.  Interferometer  Calibration  Source. 

Instrument  Response 

A  relative  instrument  response  was  obtained  each  time 

the  instrument  was  realigned  during  the  data  taking  process. 

The  calibration  alignment  process  was  repeated  about  every  2 

hours.  The  interferograms  obtained  from  the  instrument  when 

observing  the  calibration  source  were  transformed  using  the 

same  FFT  and  apodization  (discussed  later  in  this  chapter) 

routines  as  were  used  for  the  airglow  data  signal  processing 

described  above.  The  calibration  source  spectral  or 

"blackbody"  curves  were  then  averaged  together  (from  5  to  10 

frames)  to  minimize  noise  and  irregularities. 

The  tungsten  bulb  in  the  calibration  source  has  a 

blackbody  equivalent  temperature  of  2370  »K  when  operated  at 

the  specified  current  (750  mA) ,  according  to  Gilway 

Technical  Lamp  t  19821,  its  manufacturer.  The  Plank  equation 

can  be  used  to  calculate  the  spectral  sterance  L  (1)  in 

B 

watts /meter ,  for  the  assumed  blackbody  radiation  source  as, 
shown  by  Wyatt  C 19783 
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where 

h  =  6.6262  x  10-34  (Js)  (Plank’s  constant)  tWyatt  19781, 
c  =  2.9979  x  108  (m/s)  (speed  of  light)  CWyatt  19781, 

X  «=  wavelength  (m)  , 

k  =  1.3806  x  lO  23  (J/*K)  (Boltzmann’s  constant) 

CWyatt  19781, 

T  «=  absolute  temperature  (*K>. 

The  interferometer  detector  measures  energy  coming  -from 
the  interference  of  two  light  rays  of  the  same  optical 
frequency!  and  the  germanium  detector  operates  in  a  photon 
sensitive  mode.  As  a  consequence,  in  the  analysis  it  is 
desirable  to  manipulate  Plank’s  equation  into  terms  of 
wavenumber  and  quanta  (photons).  The  relationship  between 
sterance  and  photon  sterance  is 

L  =  L„X/(A<r)  .  (3.6) 

P  B 

Noting  the  relationship  between  wavenumber  and  wavelength 
and  that  Equation  3.5  is  a  density  function 

d  =  1/X  ,  (3.7) 

<fV  -  -  1/X2  flTX  .  (3.8) 

Using  Equations  3.7  and  3.8  and  converting  from  meters  to 

centimeters,  Plank’s  equation  can  now  be  shown  in  terms  of 

.  „  -2-1  -1  -1 
photons  cm  sr  sec  cm  or 
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The  photon  sterance  for  the  2370* K  tungsten  bulb 
calibration  source  Has  calculated  using  Equation  3.9  and 
corrected  -for  the  emissivity  of  tungsten  using  data  -from 
Weast  C 19771.  The  averaged  blackbody  spectrum  taken  by 
viewing  the  calibration  source,  was  then  divided  by  the 
photon  sterance.  The  resulting  curve  was  normalized  to  its 
peak  value  and  constitutes  the  relative  instrument  response 
used  in  the  rotational  temperature  calculation  model 
developed  in  Chapter  IV .  The  values  for  the  instrument 
response  curve  are  shown  as  part  of  Tables  4-1,  4-2,  4-3, 
and  4-4. 

Phase  Correction 

The  spectrum  computed  from  an  i nterf erogram  using  the 
selected  FFT  is  algebraically  complex  (contains  both  real 
and  imaginary  parts)  [Ware  19801.  The  incoming  signal 
contains  no  inherent  phase  information!  therefore,  any  phase 
angle  computed  during  the  FFT  process  is  an  artifact  of  the 
system.  A  nonzero  phase  relationship  is  caused  by  two  main 
factors.  First,  phase  shifts  occur  in  optically-compensated 
interferometers  because  the  optical  beamsplitter,  wedges, 
and  lenses  are  not  strictly  uniform  as  a  function  of 
wavenumber.  Second,  a  linear  phase  shift  occurs  during  the 
sampling  of  the  i nterf erogram  for  FFT  processing,  when  the 
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zero  path  difference  point  is  not  centered  in  the  signal 
time  window. 

The  amplitude  of  the  measured  spectrum  could  be  obtained 
by  the  magnitude  operation  (square  root  o-f  the  sum  of  the 
squares  o-f  the  real  and  imaginary  parts).  However,  this 
operation  always  yields  positive  noise  components, 
increasing  the  noise  by  f2.  This  increase  in  noise  is 
especially  harmful  when  several  -frames  o-f  data  are  signal 
averaged  (coadded)  to  help  identify  low  weak  -features  in  the 
airglow  emission  spectrum.  The  sampling  phase  shift  error 
may  also  vary  from  frame  to  frame.  Consequently,  a  method 
of  phase  correction  based  upon  the  data  within  each  frame 
must  be  used. 

The  phase  character i st i cs  of  each  data  frame  could  he 
obtained  by  Fourier  transforming  a  small  data  set  around  the 
"center"  of  the  i nterf erogram.  This  process  would  yield  a 
very  low  resolution  spectrum  from  which  the  slowly  varying 
phase  information  could  be  easily  extracted.  However,  this 
truncation  and  transform  process  is  the  same  as  convolving 
the  original  spectrum  with  the  Fourier  transform  of  the 
truncating  function.  As  described  by  Ware  119803,  this 
convolution  is  merely  a  digital  filter  operating  in  the 
f requency /phase  domain.  Hamming  t 19773  showed  that  for  a 
given  filter  width  the  minimum  noise  gain  is  obtained  when 
the  convolving  function  is  rectangular  in  shape.  The  phase 
character ist ics  are  not  related  to  the  incoming  signal  in 
any  way  so  the  rectangular  filtering  process  may  be 
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repeated  as  many  times  as  necessary  to  identify  and  remove 
the  phase  anomal ies  from  the  computed  spectrum.  When  the 
filtering  process  is  complete,  the  remaining  "real"  part  of 
the  transform  is  used  as  the  measured  spectrum  (power 
spectral  density  function). 

The  spectral  frames  were  phase  corrected  with  convolving 
filter  windows  of  3,  5,  9,  17,  and  65  sample  point  widths  in 
a  repetitive  manner.  Figure  3-3  shows  an  instrument 
uncorrected  blackbody  curve  both  before  and  after  phase 
correction,  and  Figure  3-4  shows  the  same  information  for  a 
spectral  frame.  Examination  of  Figures  3-3  and  3-4  show 
that  the  imaginary  portion  of  the  power  spectral  density 
(caused  by  chromatic  variations  within  the  instrument 
optics)  has  been  eliminated  by  the  phase  correction 
a  1  gar i thm . 


Apodization  and  Interpolation 

The  calculation  of  OH  airglow  rotational  temperatures 
was  a  major  goal  of  this  study.  The  calculation  of 
temperature  requires  the  extraction  of  relative  spectral 
line  intensities  from  the  transform  of  the  i  nterf erogram. 
The  relative  instrument  response  and  phase  correction  steps 
described  above  are  modifiers  in  signal  processing  the  line 
i n  tens i t i es . 

The  shape  of  an  individual  spectral  line  is  of  interest 
because  the  shape  will  affect  how  the  a  best  estimate  of  the 


1  i  ne 


emission  intensity  is  extracted  from  the  data.  Due  to 
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(a. )  Real 


Imag i nary 


(  b .  )  Rea 1 


Imaginary 


Figure  3-4.  Real  and  imaginary  parts  of  a  spectral  curve 
(a. )  before  and  (b. )  after  phase  correction. 


the  discrete  nature  of  the  transitions  of  the  excited  OH 
molecule,  the  spectral  lines  entering  the  interierometer  are 
discrete  in  shape  except  for  a  small  finite  width  due  to 
Doppler  and  collisional  broadening  plus  "time  windowing"  of 
each  photon  source.  The  numerical  computation  of  an  FFT 
requires  that  the  sampling  be  limited  in  length.  The 
truncation  or  multiplication  of  the  interferogram  by  a 
rectangular  "window*  in  the  time  domain  is  the  same  as 
frequency  domain  convolution  of  the  discrete  spectral 
frequencies  with  the  Fourier  transform  of  the  rectangular 
window  function.  The  Fourier  transform  of  the  rectangular 
window  is  a  (sin  x)/x  or  "sine"  function.  This 
characteristic  spectral  line  shape  is  referred  to  as  an 
"instrument  function."  The  high  side  lobe  behavior  of  the 
sine  function  (largest  side  lobe  -13  dB  down  in  amplitude) 
causes  the  various  lines  in  the  spectrum  to  interact  or  mesh 
together.  Therefore,  it  is  desirable  to  establish  a 
technique  to  suppress  the  sine  function  behavior  and  force 
the  signal  component  into  as  much  of  a  discrete  line  shape 
as  possible. 

The  process  of  suppressing  the  sine  function  side  lobes 
is  called  apodization.  The  apodization  of  the  spectral 
information  can  be  accomplished  either  in  the  i nterf erogram 
domain  by  multiplying  the  interferogram  with  another  window 
then  performing  the  FFT,  or  by  convolving  the  trans-form  of 
the  window  with  the  transformed  i nter f erogr am .  The 


convolution  of  the  spectral  data  with  the  trans-form  of  the 
chosen  window  is  the  method  used  for  this  study. 

Numerous  apodization  functions  have  been  tried  for  use 
with  Fourier  transform  spectroscopy.  In  choosing  an 
apodizing  function  the  trade-off  is  among  instrument 
function  width  ( r eso 1 ut i on ) ,  side  lobe  attenuation,  and 
computational  efficiency.  Norton  and  Beer  C19761  computed 
over  llOO  different  apodization  windows  and  plotted  each  as 
a  function  of  central  peak  width  versus  height  of  maximum 
side  lobe.  Their  study  indicated  an  optimal  boundary 
existed  between  the  two  plotted  parameters.  Figure  3-5 
shows  the  Norton  and  Beers  limit  with  some  specific 
functions  also  shown.  Vagin  119801  independently  computed 
over  3000  functions  and  also  showed  that  this  boundary 
exists.  However,  he  went  on  to  analytically  show  that  for  a 
chosen  characteristic  (either  side  lobe  suppression  or 
resolution)  an  optimal  apodization  function  can  be  computed. 
Harris  C19781  presents  an  excellent  analysis  of  many  of 
these  apodization  functions  and  shows  the  trade-offs 
associated  with  each.  Nuttall  C 19811  shows  some  corrections 
to  Harris’s  work  and  presents  additional  functions  for 
considerat ion. 

The  present  interferometer’s  full-width  half-maximum 
(FWHM)  instrument  function  is  1.8  cm  * .  The  OH  spectrum  in 
the  1  -)im  wavelength  region  has  line  separations  on  the  order 
of  10  cm  *.  This  occurrence  renders  the  prime  consideration 
in  the  choice  of  an  apodization  function  as  one  of  side  lobe 


attenuation  and  ease  of  computation, 


rather  than  maximum 


resol ut i on . 

As  can  be  seen  -from  Figure  3-5  the  Hamming  window 
(sometimes  referred  to  as  the  “minimum  2-point"  window)  is 
located  on  the  optimal  boundary.  This  -function  is 
comparatively  simple  to  compute  because  it  contains  only  two 
terms  and  provides  side  lobe  attenuation  of  -43  dB  (a  20  dB 
improvement  over  that  provided  by  the  sine  function).  The 
additional  side  lobe  attenuation  is  sufficient  because  the 
data  collected  for  this  study  have  s i gnal - t o -no i se  ratios  of 
about  lOO,  thus  placing  the  side  lobe  behavior  below  the 
noise  level.  The  Hamming  window  yields  goad  spectral 
resolution  by  providing  a  FWHM  central  lobe  of  2.7  cm 
which  is  sufficient  to  identify  the  OH  spectral 

characteristics.  This  window  also  has  an  asymptotic  side 
lobe  roll-off  of  6  dB  per  octave  [Harris  19781. 

The  Hamming  convolution  function  takes  the  form 

k  =  +  6  r 

I  =  E  I  0.53836  sincUclkiff*)) 

k— 6 

♦  0.46164  fs i nc (w [ c ( kdff  «  )+ll)+sinc<w[c(kAC«  )  -  1  3  > 1 

2 


(3.10) 


where 

«=  C  -d  «  , 

k  *  counter  for  integer  sample  numbers  centered  at 


transform  point  nearest  o», 
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—  wavenumber  at  integer  sample  in  cm  , 

=  wavenumber  at  point  where  intensity  is  being 
computed  in  cm  1 , 

=  intensity  of  transformed  data  at  integer  transform 
wavenumber  O' , 

*  intensity  of  apod  i  zed  data  at  wavenumnber  o», 

=  number  of  sample  points  per  wavenumber. 
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Figure  3-5.  Optimal  apod i zat ion  functions, 

resolution  vs.  side  lobe  attenuation  [Espy  19841. 


show i ng 


The  discrete  spectral  data  points  computed  by  the  FFT 
occurred  at  1.5-cm  intervals!  therefore,  as  can  be  seen 
from  Equation  3.10,  the  convolution  computation  was  summed 
over  t  4  sample  number*.  E<aminatlon  of  Equation  3.10  shows 
that  this  equation,  in  addition  to  suppressing  side  lobes, 
can  be  used  to  Interpolate  spectral  values  between  the 
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discrete  computed  points  of  the  FFT.  Figure  3-6  shows  the 

Hamming  instrument  -function  compared  with  the  “square 

window*  sine  -function  and  also  the  common  “triangle  window" 
2 

sine  -function.  All  three  -functions  have  been  normalized  to 
their  respective  peak  values  at  1000  -for  comparison. 


Figure  3-6.  Comparison  o-f  normalized  instrument  functions; 
sinc< - ),  sine  <-  -1,  and  Hammingl'  1  >  )  CEspy  19843. 


Line  Amplitude  Extraction 


The  i nter f erograms  were  recorded  in  “raw"  form  on  analog 
tape  during  the  measurement  campaigns.  The  analog  tapes 
wei  ater  played  back,  digitized,  formatted  by  computer  and 
stored  on  digital  tape.  The  digitized  i nterf erograms  were 
then  submitted  to  the  FFT  signal  processing  routine.  The 
power  spectral  density  information  generated  by  the  FFT  for 
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each  i n t e r f er ogr am  was  also  stored  on  digital  tape  •for  later 
pr  ocess i ng . 

The  information  of  interest  within  each  spectral  frame, 
needed  for  rotational  temperature  c a  1 cu 1  at i ons ,  was  the 
intensity  of  a  few  specific  OH  spectral  lines  (the 
temperature  model  is  explained  in  Chapter  IV).  In  order  to 
identify  the  spectral  features  of  interest,  each  of  the 
transformed  frames  of  data  was  plotted  on  paper  using  an  x-y 
recorder  and  a  simple  computer  routine  which  performed  a 
straiqht  line  connect  between  each  of  the  discrete  points 
computed  by  the  FFT.  Each  of  the  plotted  frames  was 
examined  manually,  and  the  lines  of  interest  were 
identified  by  FFT  sample  number.  The  lines  selected  were 
bright  lines  within  each  OH  band  which  were  far  enough  apart 
as  to  not  be  contaminated  by  other  features  and  situated  at 
wavelengths  of  minimal  atmospheric  absorption  by  H^O  and 
CO^.  Examination  of  the  plots  showed  that  the  position  of 
each  line  never  varied  from  frame  to  frame  more  than  ±1  FFT 
sample  number.  The  selected  FFT  sample  point  was  then 
recorded  as  the  initial  position  of  each  of  the  spectral 
lines  of  interest. 

The  extraction  of  spectral  line  intensities  was 
accomplished  using  the  Hamming  apodization  function  to 
interpolate  between  the  discrete  FFT  data  points.  Allowing 
far  additional  error  in  the  chosen  position  of  each  line, 
the  Interpolation  routine  was  operated  in  a  sample  number 
window  13  data  points  around  that  manually  chosen  line 


6  1 


position.  The  line  intensity  was  then  computed  at  0.01 
sample  number  increments  <0.013  cm  *  increments)  over  the 
data  point  window  using  the  Hamming  function.  The  routine 
saved  the  maximum  value  found  within  the  sample  window  as 
the  line  amplitude  for  that  spectral  feature. 

The  instrument  calibration  function  and  the  correct 
spectral  line  amplitudes  could  now  be  used  in  the 
temperature  model  tor  extraction  of  OH  rotational 
temperatures  as  a  function  of  each  i nterterometer  scan. 
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CHAPTER  IV 

OH  ROTATIONAL  TEMPERATURE  MODELING 

Introduct ion 

Hydroxyl  is  a  minor  atmospheric  constituent,  residing 
in  a  thin  layer  (about  7  km  thick)  at  an  altitude  near  87  km 
[  Baker  et  al.  19853.  Although  minor  in  concentrat i on 
(between  10^  and  106  mo  1  ecu  1 es / cm3  at  night,  Baker  [19783) 
it  is  the  major  atmospheric  near-infrared  airglow  radiator 
at  night.  The  radiation  generated  by  OH  occurs  in  spectral 
bands  known  as  the  Meinel  bands,  since  their  discovery  and 
identification  by  Meinel  [19503.  After  their  discovery,  the 
measurement  of  these  bands  has  been  of  great  interest  to  the 
atmospheric  science  community.  OH  airglow  emission  band 
measurements  contain  in-formation  about  mesospheric 

populations  and  temperatures ,  and  as  such  provide  insight 
into  the  dynamics  of  the  entire  middle  atmospheric  region. 
A  rotat ional -temperature  model  is  developed  in  this  chapter 
similar  to  the  technique  used  by  Hill  et  al .  [19793.  The 

field  measurements  are  then  fit  to  the  model  in  a  least- 
squares  sense.  The  quality  of  the  fit  will  also  determine 
error  bounds  on  both  relative  band  intensity  and  absolute 
rotational  temperature. 
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The  Excited  OH  Radical 

The  hydroxyl  radical  is  a  diatomic  heteronuclear 
molecule  and  has  a  nonzero  electric  dipole  moment.  The 
molecule,  therefore,  has  the  ability  to  be  easily  excited 
and  radiate  electromagnetic  radiation.  The  OH  Meinel 
radiation  bands  in  the  infrared  are  caused  by  vibration- 
rotation  transitions  within  the  ground  electronic  state  of 
the  OH  molecule.  (The  electronic  absorption  and  emission 
spectra  o-f  OH  occur  in  the  ultraviolet.) 

Hydration  o-f  ozone  and  perhydroxyl  reduction  are  the 
primary  processes  -for  the  creation  of  vibrational  ly-exc  i  ted 
OH  in  the  earth’s  upper  atmosphere  [Baker  19781.  The 
v i brat i onal 1 y-exc i ted  states  of  OH  are  quantized  and  are 
populated  according  to  certain  dipole  selection  rules 
described  by  Hertzberg  [19713.  The  quantum  numbers 
associated  with  the  vibrational ly-excited  states  take  on 
values  v  *  0 , 1 , 2 , 3 , . . . , 9  and  the  allowed  transitional 
changes  during  emission  are 

v* -v* •  =  6V  =  1,2,3,  . . . , 9  ,  (4.1) 

where  v'  is  the  upper-state  energy  level  and  v’ *  is  the 
lower-state  energy  level.  The  Meinel  radiation  bands  are 
known  by  the  value  of  6v  associated  with  the  transition. 
The  band  sequences  measured  in  this  study  are  Av=2  (4,2  and 
3,1)  and  iv«3  (8,3  and  7,4)  transitions  since  they  have 
energy  differences  in  the  near - t nf r ared  region.  A  typical 


measured  spectrum  is  shown  in  Figure  4-1. 


Figure  4-1.  Measured  OH  spectra  showing  (4,2),  (3,1), 

(8,5),  and  (7,4)  Meinel  bands. 

Each  band  sequence  shown  in  Figure  4-1  is  a  complex 
group  of  spectral  lines.  The  complexity  is  due  to 
transitions  in  molecular  rotational  energy  in  addition  to 
the  changes  in  vibrational  levels.  These  angular  momentum 
changes  modulate  the  vibrational  transitions  causing  the 
intricate  spectral  structure.  Appendix  A  includes 

additional  information  on  OH  radical  transition  theory. 

Rotational  Temperature  Model 


In  extracting  temperatures  from  OH  spectral  data  it  is 
assumed  that  the  OH  molecular  population  is  in  rotational- 
state  thermal  equilibrium.  The  population  distribution  over 
the  different  quantum  numbers  can  then  be  described  by  the 


Maxwe 1 1 -Bo  1 tzmann  distrioution  taw  (Baker  1978], 


N  .  ..  =  number  of  molecules  in  vibrational  level  v’  and 

v'J* 

rotational  level  J*  at  temperature  T, 

w  f  =  2J’>1  or  the  statistical  weight  state  J '  with 

its  <  2J '  ♦  1 )  -f  o  1  d  degeneracy ,  . 

C  ■  a  constant.  equal  to  N  ./Q  where  Q  is  the 

v'  v'  r  r 

partition  •function  for  rotational  level  v, 

E  t  *  energy  of  the  rotational  state,  and  equal  to 

F  .  (J*  1  he.  where  F  .  <J’  )  is  the  term  value  for 
v*  v’ 

the  upper-rotat ional  state, 
k  *=  Boltzmann’s  constant. 

- 1  -3 

The  volume  emission  rate  in  photons  sec  cm  of  a 
central  line  is  given  by  Mies  11974]  as 


where 

I  ...  .  *  intensity  of  emission  from  upper-state 

v'J'  to  lower-state  v,,J'', 

A  ...  .....  “  Einstein  transition  probability  of 

v'J'  ,  v’  *  J '  ' 

spontaneous  emission  from  state  v'J'  to 
state  v ’ ' J  * ’  . 

Inserting  Equation  4.2  into  Equation  4.3  gives  the  volume 
emission  rate  of  each  rotational  line  arising  from  a 
transl t ion  from  an  upper-state  v'J’  to  a  I ower -state  v’  '  J  '  '  . 


<4. 4) 


|  The  computation  of  a  rotational  temperature  from 

Equation  4.4  requires  a  value  -for  the  absolute  line 
intensity  Iv,j,  however,  the  ratio  of  two  relative 

line  intensities  can  be  taken  and  a  temperature  computed 
[Ware  19801  without  the  need  -for  precise  absolute  spectral 
radiance  values.  However,  the  ratio  technique  proves  to  be 
a  nonlinear  (logarithmic)  process  and  when  several  line  pair 

i 

| 

I  temperatures  within  a  band  are  calculated,  the  averaging  o-f 


i 


l 

! 

i 

I 

! 


1 


the  pair  temperatures  to  obtain  a  true  band  temperature  is 
difficult.  The  line  intensities  are  extracted  -from  data 
whose  noise  characteristics  are  spectrally  -flat  (white 
noise)  and  gaussian  distributed  [Ware  19803.  The  di-fficulty 
comes  about  in  the  nonlinear  temperature  computation  because 
the  uncertainty  o-f  the  calculation  is  not  simply  a  linear 
extension  o-f  the  input  data  uncertainty.  The  combining  o-f 
the  line  pair  temper atures ,  therefore,  cannot  be  computed 
using  traditional  averaging  techniques. 


A  model  of  the  Boltzmann  distributed  data  was  developed 
to  eliminate  the  problems  described  above.  The  model  would 
describe  all  the  lines  that  could  occur  within  a  band  and 
then  be  fitted  to  the  measured  data  using  a  1  east -squares 
technique.  The  mathematical  model  is 


M  (A  , B  ) 
i  *  t 


A  R  I  exp  (  -F  B  ) 
»  i  1 _ i  • 


n  C I  exp ( -F  B  ) 3 
i  j  • 


(4.3) 
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67 


where 


M.  =  the  value  of  the  model  at  wavenumber  i, 

i 

A  *  relative  total  integrated  band  intensity, 

• 

3 

I  «  A(2J'+1)  =  cr  S,  which  is  the  Einstein  coefficient 

l 

times  the  J*  state  degeneracy  or  the  wavenumber 
cubed  times  the  line  strength, 

B  =  hcl j*T,  where  T  is  rotational  temperature  in  *K, 

*  term  value  for  upper-state  transition, 

R.  =  relative  instrument  response  at  wavenumber  i, 

E  m  summation  over  all  lines  within  the  band,  is  a 
scaling  factor  to  ensure  the  entire  band  intensity 
is  contained  within  the  A  term. 


The  model  shown  above  is  similar  to  the  intensity 
calculation  shown  in  Equation  4.5.  The  difference  is  that 
the  model  pertains  to  an  entire  band  of  transitions  and  is 
modified  by  R.  to  appear  like  the  measured  data.  All 
quantities  within  the  model  are  now  known  except  the  band 

intensity  A  ,  and  the  temperature  term  B  .  Within  each  OH 

•  a 

band,  the  measured  data  points  can  be  subtracted  from  the 
model  at  the  same  wavenumber  and  any  differences  will  be 

contained  within  small  changes  of  A  and  B  times  the 

•  a 

derivative  of  the  model  at  points  A  and  B  . 


a 


a 


D,  -  M,  -  a A 
i  i  a 


3M  (A  ,  B  > 

i  a  a 

a  a 


]  ♦  AB.[ 


3  M  (A  ,  B  ) 

i  a  a 

a  b 


]  . 


(4.6) 


a 


a 


where  is  the  measured  data  point.  Equation  4.6  can  be 

shown  for  <j)  data  points  within  a  band  in  matrix  form  as 


V’1 


./I 


$ 

I 


$ 

4*l! 


$ 

$ 

ijft 


,Vl 

•  7*1 


<D  -  M  ) 

l  t 


<D  -  M  ) 
j  j 


( - ->  ( - -) 

3A  3B 

«  o 


3M  3M . 

( - — )  ( - 

3  A  3  B 

o  « 


(•4.7) 


where  the  left  matrix  is  the  difference  between  each 
measured  line  intensity  within  a  band  and  the  line  intensity 
as  predicted  by  the  model.  The  right  hand  side  of  Equation 
4.7  must  be  equal  to  the  preceding  difference. 

The  computation  o-f  Equation  4.7  requires  the  solution 
for  each  of  the  derivatives 


R.  »  exp  (  -F  .  B  ) 
lilt 


(4.8) 


and  , 


A.  R  I .  exp ( -F  ,  B  ) 
l  i  i  i  • 


f  IE  (IF  exp { -F  B  > ) 1 
.  j  j  j  • 

J 


"  N  Fj  . 


(4.9) 


where  N  =  £  I  exp(-FjB  )  or  the  normalization  factor.  All 
J  j  • 

the  terms  in  Equation  4.7  are  computed  and  submitted  to  a 
1  east -squares  fitting  routine,  developed  by  Lawson  and 

Hanson  C 19743  which  returns  the  values  of  4 A  and  aB  .  The 

•  • 

least-squares  routine  is  used  to  solve  problems  of  the  form 


Ax  -  B 


(4.10) 


*v 


where 


A  =  the  derivative  matrix  of  Equation  4.7, 
x  =  the  unknown  or  A  matrix, 

B  “  the  difference  matrix. 

To  begin  the  process  a  starting  temperature  was  chosen 
and  all  terms  computed.  The  result  of  this  first  fit  is  to 

determine  an  appropriate  initial  value  for  A  (A  is  merely 

•  » 

a  relative  intensity  term).  Once  the  initial  values  for  B 

» 

and  A  were  ascertained,  the  least-square  fitting  routine 
was  operated  in  a  loop,  the  derivative  matrix  vectors  were 
each  scaled  to  unity  value  before  subm ission  to  insure 
convergence  of  the  fit,  and  each  time  the  routine  would 

return  values  for  AA  and  aB  .  The  A  values  were  added  to 

a  • 

the  previous  values  for  A  and  B  ,  the  model  values 

•  • 

recomputed  and  resubmitted  to  the  routine.  The  value  of  aB 

• 

was  tested  after  each  iteration  until  the  temperature  change 
(temperature  is  part  of  B  )  from  the  previous  fit  was  less 
than  0.5*  K.  When  the  1  east -squares  routine  converged  to 

the  0.5  *K  limit  the  values  for  A  (relative  band  intensity) 

§ 

and  T  (the  temperature  part  of  B  )  were  stored  along  with 

t 

the  time  and  date  of  the  frame. 

The  constant  terms  used  in  the  model  were  computed  by 
Espy  C1984J  using  the  molecular  constants  measured  by  Covan 
and  Foster  119821  and  OH  dipole  moments  derived  by  Werner  et 
al.  119831.  The  values  of  these  terms  are  given  in  Tables 
4-1  through  4-4.  The  tables  show  the  wavenumber,  a  typical 
relative  instrument  response,  the  term  value  for  the  upper  - 
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TABLE  4- 

L  i  ne 

1.  Molecular 

Wavenumber 

( cm- *  ) 

data  -for 

Re  1  at i ve 
Response 

OH  4 v=2 ,  (4, 

Term  Value 

(cm- * ) 

2)  band. 

Line  Strength 

Pi 

(4)  * 

6159-634 

0.  147 

13421 . 921 

O. 23415100E+13 

Pi 

(4)  * 

6174-994 

0.166 

13514. 125 

0. 19663726E+13 

Pi 

<3>* 

6200.353 

0.  199 

13320. 758 

O. 1 6226 1 62E+ 1 3 

Pt 

(3)  * 

6219.070 

0.224 

13428.966 

0. 1 322644 1E+ 13 

Pi 

( 2 )  •#■ 

6238. 114 

0. 247 

13248.923 

O. 88080260E+12 

Pi 

(  2)  * 

6260.766 

0.275 

13377.615 

O. 70291000E+12 

Qt 

<  3 ) 

6301.517 

0.  326 

13421 . 922 

O. 25067365E+12 

Qi 

(2) 

6309.948 

0.336 

13320.758 

O. 71495030E+12 

Qt 

<  1  ) 

6315.828 

0.344 

13248.923 

O. 1 16S1950E+ 13 

Rt 

(  1  ) 

6368.381 

0.411 

13428.966 

O. 60586570E+12 

R  i 

(  1  ) 

6387.661 

0.440 

13320.758 

O. 74062960E+12 

R  i 

(2) 

6411.112 

0.473 

13421 . 921 

0. 12600000E+ 13 

<* 

1  i  nes 

used  in  -fit) 

TABLE  4 

-2  .  Mo  1  ecu  1 ar 

data  for 

OH  iv=2,  (3, 

1 )  band . 

L  i  ne 

Wavenumber 

Re  1  at i ve 

Term  Value 

Line  Strength 

<  cm- 1  ) 

Response 

( cm- 1  ) 

Pi  (4)* 

6480. 230 

0.572 

10352.446 

O. 27218330E+ 13 

P.  (4)* 

6495. 579 

0.593 

10443.293 

O. 22819520E+13 

Pi  <  3  )  * 

6522. 332 

0. 626 

10247.071 

0. 18861 154E+ 13 

Pi  (3) * 

6541 . 239 

0. 650 

10354.206 

O.  15350037E+  1  3 

Pi  ( 2) * 

656 1 . 409 

0. 679 

10172.300 

0. 1 02406 17E+ 13 

Pt  ( 2  )  * 

6584 . 562 

0.712 

10300.410 

O. 81566820E+12 

Qt  (3) 

6627.706 

0. 770 

10352.446 

0. 57940270E* 12 

Qt  (2) 

6636. 171 

0. 782 

10247.071 

O . 82784900E+ 1 2 

Qt  (  1  ) 

6642. 081 

0. 789 

10172. 300 

0.  1 355684 4E+  13 

Rt  (  1  > 

6697. 063 

0.867 

10354.206 

0. 35196955E+12 

R  i  <  1  ) 

6716.853 

0. 892 

10247.071 

0 . 86348040E  *  1 2 

R i  (2) 

6741 . 555 

0.931 

10352.446 

O. 1470351 1E+13 

( *  lines  used  in  -fit) 
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TABLE  4-3.  Molecular  data  -for  OH  &v  =  3,  (8,5)  band. 


L  i  ne 

Wavenumber 

<  cm-1  > 

Relative 

Response 

Term  Value 

(cm -  *  ) 

L i ne  Strength 

Pi 

<5) 

7355. 155 

1 .  OOO 

24163.096 

O. 54607000E ♦  1  3 

Pt 

<  5 ) 

7569.641 

0. 997 

24248. 912 

0. 4  719471 OE ♦ 1 3 

P« 

(4) 

7600. 862 

0.992 

24055.028 

0.419139?0Ef  13 

Pt 

(4) 

7617.468 

0. 990 

24152. 797 

O. 35340530E ♦  1  3 

Pi 

(3)  * 

7642.343 

0. 984 

23971 . 237 

0. 2920 1 550E  ♦  1  3 

Pt 

(3)  * 

7661 . 252 

0. 983 

24083 . 736 

O.  23944700E  ♦  1  3 

Pi 

(2)  * 

7679.631 

0. 978 

2391 1 . 382 

O.  1 3952644E  ♦  1  3 

Pt 

(2)* 

7700.808 

0. 972 

24042. 1 13 

0.  1 2846050E ♦  1  3 

Qi 

<  3 ) 

7726. 136 

O.  960 

24053. 028 

O. 95265240E  ♦  12 

Qt 

(2) 

7739.286 

0.932 

23971 . 237 

0.  1 3444 1 94E  ♦  1  3 

Qi 

(  1  ) 

7748.311 

0.945 

23911 . 582 

0. 2 1 738030E ♦ 1 3 

Rt 

(  1  ) 

7791 . 408 

0.900 

24803 . 736 

O.  1  146  38  2 3E  ♦  1  3 

Ri 

<  1  ) 

7808. 166 

0.874 

23971 . 237 

0.  1 39004  38E ♦  1  3 

Rt 

(2) 

7823.076 

0.857 

24053 . 028 

O.  23699920P  ♦  1  3 

<* 

1 1 nes  used  In  fit) 

TABLE  4-4.  Molecular  data  for  Av=3,  (7,4)  band. 


L  i  ne 

Wavenumber 

( cm-‘  ) 

Re  1  at  1 ve 
Response 

Term  Va 1 ue 

<  cm- •  ) 

Line  Strength 

Pi 

(5)* 

8049.708 

0.786 

21763.099 

0.66072740E+13 

P« 

(5)  * 

8064. 112 

0.  779 

21847. 340 

O.  56979800E+1 3 

P. 

(4  >* 

8096.397 

0.765 

21649. 131 

0. 40669700E+  13 

P» 

(4)* 

8113. 112 

0. 754 

21745. 519 

0. 42663300E*  1  3 

P. 

(3)* 

8138.903 

0. 740 

21560. 824 

0. 3530245OE ♦ 1 3 

Pt 

(3)  * 

8158.204 

0.  725 

21672. 328 

0. 28898620E+  1  3 

Pt 

<2>  # 

8177.237 

0.711 

21497.994 

0. 1 9277462E ♦ 13 

Pt 

<  2 )  * 

8199.208 

0. 696 

21628. 204 

O. 15491604E+ 13 

Qi 

<  3 ) 

8227. 210 

0. 680 

21649. 131 

O. 1 1407934E* 13 

Qi 

(2)* 

8240.067 

O.  672 

21560.824 

0.  16141494E  +  13 

Qt 

<  1  )  * 

8249.071 

0 . 668 

21497.994 

O. 26171900E+  13 

Rt 

<  1  )* 

8294 .714 

0. 646 

21672. 328 

0.  1 3787860E ♦  1  3 

Ri 

(  1  )  * 

8311. 901 

0. 642 

21560.824 

0.  1 676 1 780E ♦  13 

Ri 

(2)* 

8328.374 

0 . 635 

21649. 131 

0. 28788900E  +  13 

<#  lines  used  In  fit) 


state  transition,  and  the  theoretical  line  strength  tor  each 
1  ire.  The  aa,  emimher  shown  is  the  average  ot  the  two  A- 
doubled  lines;  the  term  value  is  also  the  average  o t  the  two 
lines,  and  the  line  strength  given  is  the  sum  at  the  two 
line  strengths.  The  sum  o+  the  line  strengths  was  used 
because  the  measured  line  is  actually  the  sum  o-f  the  A- 
doubled  pair.  Data  tor  all  the  strong  lines  are  listed  tor 
use  in  the  nor ma 1  i z a t  1  on  tactor  ot  the  model  but  only  those 
with  were  used  in  the  tit.  The  reason  all  were  not  used 
in  the  tit  was  that  many  were  determined  to  have  been 
contaminated  b  water  vapor  absorption  or  by  other  emission 
1  i  n e s  C  R o  .■  c  h ou  r  d  h u  r  /  19831. 

Error  Analysis  and  Testing  ot  Model 

The  model  used  here  is  over-determined  with  only  two 
unknowns  and  at  least  tour  equations  (in  (8,5)  band  tit, 
more  in  other  bands)  and  as  such,  the  additional  intormation 
can  be  used  to  estimate  the  accuracy  ot  the  1  east -squares 
tit.  Examination  ot  Equations  4.7  and  4.10  show  that  when 
the  titting  routine  is  complete,  the  derivative  matrix  A  is 
lett  as  a  2X2  upper - tr i angu 1 ated  matrix.  The  rest  ot  the 
terms  in  the  A  matrix  are  ot  little  signiticance  because  the 
1  east -squar es  routine  has  manipulated  the  data  into  the 
upper-triangulated  portion  in  order  to  solve  the  equations. 
The  upper- triangulated  matrix  now  takes  on  the  torm 


' 

i 


# 


Z3E 


Mhart 


H  •  upper - t r t an^u 1 ated  coefficient*  returned  in  the 
derivative  matrix  A. 

Upon  conpletion  of  the  -fitting  routine,  the  values  left 
in  the  matrix  x  must  be  equal  to  zero,  except  for  noise  in 
the  fit,  and  the  upper  two  values  in  the  vector  B  must  also 
be  equal  to  zero,  also  except  for  noise,  because  the  system 
of  equations  has  been  solved.  But  because  noise  does  exist 
in  the  fit  then  some  residual  remains.  Therefore,  the 
equation  system  is  left  in  the  form 

ft  x  -  n  ,  (4.12) 

where  n  is  the  noise  vector. 

In  order  to  calculate  the  accuracy  of  the  models’  fit  to 
the  two  system  unknowns  (temperature  and  band  Intensity)  it 
is  necessary  to  obtain  the  covariance  matrix  which  contains 
the  variance  of  each  parameter  in  the  1  east -squares  fit.  It 
can  now  be  shown  by  manipulation  of  Equation  4.12  that  the 
covariance  matrix  is  (Lawson  and  Hanson  19741 

<x  xT>  ■  R  <n  n^>  RT 


(4. 13) 
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whert 

<  x  xT>  ”  covarlanci  matrix, 

R  *  Inverse  of  R, 

T  indicates  the  transpose  operation, 

<n  nT>  *  the  variance  of  white  uncorrelated  noise. 

The  noise  in  the  system  fit  is  assumed  to  be  white  and 

uncorrelated  because  all  the  mathematical  manipulations 

performed  were  linear  in  nature. 

A  measure  of  the  noise  in  the  system  is  returned  by  the 

1  east -squares  fitting  routine.  The  returned  parameter  is 

the  Euclidean  norm  of  the  residual  vector  and  is  called 

R  .  The  variance  of  the  system  fit  can  be  computed  from 

norm 

R  as 

norm 

d2  -  (R  2) / (M-K)  ,  (4. 14) 

norm 

where 

2 

d  *  variance  of  fit, 

M  *  number  of  equations  in  the  fit, 

K  ■  number  of  unknowns  in  the  fit. 

The  covariance  matrix  can  be  solved  in  terms  of  the  variance 
of  the  model's  fit 

<x  xT>  =  d2  I  <RT  R)  ,  (4.15) 

where  Z  is  the  identity  matrix. 

The  standard  deviation  for  each  of  the  fit  variables  can 

be  computed  from  Equation  4.15  in  terms  of  the  original 

2 

variables  of  the  residual  matrix  R  and  the  variance  d  . 
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2x2- 


( r 


1.1 


where  C  is  the  standard 
A 

function. 

The  standard  deviation  on 
part  of  B)  is 


♦  r  ) 

- ^ -  (ff  2) 


r2,2> 


deviation  of 


the  variable 


B 


2 


(4.16) 


the  intensity 


(temperature  is 


Referring  to  Equation  4.5  for  the  relationship  between 
temperature  T  and  the  variable  Bt  the  standard  deviation  of 
the  temperature  can  be  calculated 

dT  **  (ffgAc*)  /  (B#2)  ,  (4.18) 

where  is  the  standard  deviation  of  the  temperature  as 
calculated  by  the  model. 

The  model  presented  herein  provides  a  method  for 
calculating  relative  band  intensity,  rotational  temperature, 
and  the  standard  deviation  on  both  parameters.  Appendix  D 
Presents  the  computer  programs  which  implement  the  data 
processing  algorithms  developed. 

As  verification  of  the  quality  of  the  model,  the  process 
outlined  was  used  to  compute  the  Intensity  and  temperature 
of  synthetic  spectra  generated  by  Espy  t 19843.  The  signal- 
to-noise  ratio  of  the  synthetic  spectra  was  varied  from  a 
low  of  2  to  a  high  of  20  and  then  each  was  submitted  to  the 


Figure  4-2.  Model  temperature  teat  results  on  synthetic 
spectra  vs.  signal -to-noise  ratio  with  standard  deviation 
shown  as  bars  CEspy  19843. 

model  for  processing.  Figures  4-2  and  4-3  show  the  results 
of  the  test.  As  can  be  seen,  the  mean  Intensity  and 
temperature  of  the  model  closely  track  the  actual  parameters 
even  in  very  hlgh-nolse  environments.  The  error  bars  shown 
on  both  figures  indicate  a  1  sigma  uncertainty  associated 
with  the  calculation,  and  as  the  noise  increased,  the 
uncertainty  of  the  calculation  increased  as  expected.  The 
tests  were  conducted  on  a  single  frame  of  data!  therefore, 
an  improvement  could  be  made  by  averaging  frames  at  the  cost 
of  temporal  resolution  degradation.  The  signal -to-noi se 


ratio  of  the  measured  data  was  usually  about  100, 


which  is 
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CHAPTER  V 
RESULTS 


Intraduct  ion 


The  high-throughput,  narrow  field  of  view 
i nterf erometer -spec trometer  was  taken  to  several  airglow 
observation  sites  in  the  western  United  States  in  the  spring 
of  1983.  Over  several  months  of  observation,  one  of  the 
brightest  airglow  structure  events  was  seen  and  recorded  on 
June  15,  1983,  thus  fulfilling  the  primary  goal  of  this 
study.  The  observation  was  made  from  Sacramento  Peak,  New 
Mexico.  The  site  is  located  at  105,48’16“  west  longitude, 
32*47'57"  north  latitude  at  an  elevation  of  9570  feet.  The 
results  presented  in  this  chapter  were  taken  from  this  site 
between  June  13  and  June  15,  1983.  A  complete  catalog  of 
the  results  computed  from  the  1 nterf erometer  data  is 
presented  in  Appendix  C. 

Background 


The  1 nter f erometer  with  an 
isocon  camera  mounted  on 
i n ter f erometer  telescope  < see 
collect  data  at  the  observat 
camera  was  mounted  on  its  own  tr 
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m ore  bulky  lnttrf*rom*tfr.  The  infrared  cameras  were 
supplied  by  the  University  of  Southampton,  England.  Taylor 
(1983-841  operated  the  camera  equipment  and  provided 
experience  and  expertise  in  the  search  for  airglow  structure 
events.  The  cameras  were  the  “eyes*  for  the  interferometer 
in  determining  if  any  OH  airglow  structure  was  present. 
Once  an  area  of  OH  airglow  structure  was  located  in  the  sky 
the  earner a- i nter f erometer  system  was  positioned  to  view  that 
area.  The  data  from  the  two  coaligned  systems  were  used  to 
correlate  the  viewed  and  calculated  intensities  and 
temperatures . 

An  infrared  radiometer  CHuppi  19761  was  also  used  to 
gather  trend  and  absolute  intensity  data  of  the  OH  activity. 
The  radiometer  looked  in  the  zenith  and  has  a  field  of  view 
of  9*.  The  radiometer  spectral  bandpass  was  centered  at 
1.93  Mm  wavelength. 

Infrared  Isocon  Camera  Results 

The  Isocon  camera  systems’  response  is  in  the  near- 
infrared  from  about  700  to  890  nm.  The  response  curve  and 
the  OH  transitions  which  occur  within  this  bandpass  are 
shown  in  Figure  1-3. 

Figure  9-1  shows  a  video  frame  taken  with  the  large, 
independent ly-mounted  isocon  camera  at  8: 19  hrs.  UT  on  June 
19,  1983.  The  field  of  view  for  the  photo  is  28*  vertical 
and  37*  horizontal  at  an  azimuth  of  about  320*  and  the 
bottom  of  the  picture  beginning  at  10*  elevation. 
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Examination  of  this  -frame  shows  7  distinct  bright  OH 
emission  bands.  The  bands  are  moving  northward  in  a 
direction  normal  to  the  bands  ( mov i ng  towards  1 ower  right 
hand  corner  of  frame).  These  structural  bands  extended 
across  the  entire  northern  half  of  the  sky.  The  structure 
depicted  in  Figure  5-1  was  observed  and  recorded  beginning 
at  7:30  hrs.  UT  continuing  until  10: 15  hrs.  UT  on  this  date. 

Using  the  video  record,  Taylor  C1983-843  calculated  the 
apparent  wavelength  and  velocity  of  the  band  structure,  with 
techniques  described  by  Hapgood  and  Taylor  t 19823.  Using 
the  location  and  altitude  of  the  observation  site,  radius  of 
the  earth,  viewing  angles  of  the  camera,  and  assuming  that 
the  band  structure  res'ded  on  a  spherical  shell  at  an 
constant  altitude  of  87  km,  the  apparent  temporal  wavelength 
derived  was  24±  1  km,  and  the  apparent  period  was  14±1 
minutes.  Based  on  these  two  values  the  apparent  velocity  is 
28i 2  meters /second  which  agrees  closely  with  other  similar 
observations  CTaylor  et  al .  19803. 

The  small  isocon  camera  which  was  mounted  on  and 
coal igned  with  the  i nterf erometer ,  was  operated  during  the 
same  time  frame  as  the  large  camera.  The  field  of  view  of 
the  small  camera  is  nearly  square  being  13*  vertical  by  15* 
hor i zontal . 

Figures  5-2  through  5-4  show  a  sequence  of  video  frames 
beginning  at  7:32  hrs.  UT  and  taken  at  8  minute  intervals 
with  the  small  camera.  Each  frame  again  shows,  in  more 
detail,  the  OH  emission  structure.  The  "X"  in  each  frame 


Figure  3-2.  Isocan  IR  camera  photo,  day  166,  7:32  hrs.  UT , 
■field  o-f  view  »  13*  vertical  13*  horizontal,  1  nterf erometer 
looking  at  "X*  on  a  bright  band. 


Figure  3-3.  Isocan  IR  camera  photo,  day  166,  7:40  hr*.  UT , 
field  of  view  -  13*  vertical  13*  horizontal,  Interferometer 
looking  at  "X*  on  a  dark  band. 
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marks  the  center  of  the  interferometer’s  field  of  viev 


within  the  video  -frame.  The  interferometer's  field  of  view 


is  0.8*  full  field  and,  therefore,  is  about  0.4  inches  in 


diameter  on  each  frame.  This  small  field  of  view  permits 


the  interferometer  to  resolve  the  bright  and  dark  portions 


of  the  structure.  The  i nterf erometer  is  viewing  at  an 


elevation  angle  of  17*  and  an  azimuth  angle  of  328*  in  the 


sequence.  The  three  photos  show  one  of  the  periods  where 


the  i nterf erometer  was  able  to  view  a  bright,  a  dark,  and  a 


bright  band  in  sequence  as  the  structure  moved  through  the 


field  of  view. 


Figures  5-5  and  5-6  show  another  similar  sequence  of 


video  frames  observed  later  in  the  night  beginning  at  8:31 


hrs.  UT  at  an  interval  of  11  minutes.  This  later  series  of 


frames  again  shows  the  interferometer  viewing  a  dark  then  a 


bright  emission  band  as  the  structure  moves. 


i nterf erometer  was  viewing  at  an  elevation  angle  of  15.5* 


and  an  azimuth  angle  of  340*  in  these  last  two  figures. 


Taylor  11983-843  compared  the  bright  and  dark  bands 


recorded  by  the  small  isocon  camera  with  the  OH  (3,1)  band 


intensity  plots  computed  using  the  1 nter f erometer  data 


(shown  later  in  Chapter  V  and  in  Appendix  C).  In  the  time 


period  from  7:30  to  10:15  hrs.  UT  the  two  independent 


records  show  that  for  each  bright  or  dark  band  depicted  in 


the  video  data  a  corresponding  increase  or  decrease  in 


intensity  is  also  shown  in  the  i nterf erometer  data.  This 
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comparison  showed  16  distinct  points  of  correlation  between 
the  video  -frames  and  the  interferometer  intensity  plots. 

Infrared  Radiometer  Data 

The  radiometer  was  used  to  show  trend  and  absolute 
radiance  values  for  the  OH  infrared  activity  as  viewed 
within  a  9*  field  of  view  in  the  zenith.  Radiometer  data 
are  included  to  show  for  the  night  before  (June  14,  1983) 

and  the  night  of  (June  15,  1983)  the  recorded  structure,  OH 
airglow  emission  intensity  exhibited  some  unusual  trends  and 
modulations.  Figure  5-7  shows  the  radiometer  data  for  three 
specific  days  and  is  calibrated  in  k i 1 o-Ray 1 e ighs  (kR)  of 
intensity  of  OH  (3,1)  band  emission.  The  scale  has  been 
shifted  to  show  the  three  days  all  on  one  chart. 

Day  162  is  a  typical  curve  of  near - i nf r ar ed  OH  activity 
and  was  taken  from  White  Sands  Missile  Range  (4300  feet 
elevation)  and  located  about  30  miles  from  the  Sacramento 
Peak  site.  As  can  be  seen  from  the  Day  162  curve,  the 
intensity  of  the  1.53  Mm  radiation  steadily  decreases  after 
local  sunset  (sunset  occurred  about  3:30  hrs.  UT )  and 
throughout  the  night. 

The  Day  165  radiometer  curve,  however,  shows  a  markedly 


d if f erent 

trend 

i  n 

the 

OH  activity.  The 

pos  t -sunset 

decrease 

starts , 

but 

at 

about  two  hours 

before  local 

midnight 

(or  5:00 

hrs . 

UT) 

the  OH  activity 

dr amat i cal ly 

i ncr eased 

.  After 

the 

peak  , 

the  OH  activity  remains  high  but 

also  shows  some  modulation  that  could  be  interpreted  as 


distinct  large-scale  wave  structure  passing  overhead.  The 
period  of  these  intensity  waves  ranges  -from  30  to  60  minutes 
with  a  modulation  (or  contrast  ratio)  of  from  3  to  10 
percent.  The  waves  were  not  distinguishable  in  the  zenith 
with  the  infrared  camera  equipment  because  the  zenith 
emission  intensity  is  too  faint  and  their  large  size 
exceeded  the  camera’s  field  of  view.  A  camera  was  lowered 
in  elevation  in  a  attempt  to  view  the  waves  (lower 
elevations  view  the  structure  obliquely  thus  the  waves 
appear  compressed  together  and  brighter  due  to  the  van  Rh i j n 
effect)  but  no  structures  were  d i st i ngu i shab 1 e  by  the  camera 
throughout  the  entire  night. 

The  radiometer  curve  in  Figure  5-7  for  day  166  also 
shows  enhanced  OH  activity.  Although  not  as  intense  as  the 
previous  night,  a  similar  pre-midnight  maximum  is  observed. 
After  the  pre-midnight  maximum,  the  decrease  is  slower  than 
normal  and  again  some  distinct  wave  structure  is  observed. 
The  structure  has  periods  of  about  25  minutes  and 
modulations  in  intensity  of  about  5  percent.  Shortly  after 
local  midnight  (7:00  hrs.  UT ) ,  the  moon  had  set  far  enough 
to  lower  the  sensitive  camera  to  the  horizon  to  search  for 
observable  structure.  The  measurements  gathered  throughout 
the  remainder  of  the  night  with  the  cameras  and 
i nter f erometer  also  exhibit  high-contrast  OH  emission 


struc  ture 


Direct  Comparison  of  Interferometer 


and  Radiometer  Intensities 

Prior  to  moonset  on  UT  days  165  and  166  the 
interferometer  was  pointed  to  the  zenith  and  was  viewing  the 
sky  coincident  with  the  radiometer.  The  re-fore  a  direct 
c ompar ison  is  made  between  the  two  records.  Figure  5-8 
shows  the  OH  <3,11  band  intensity  -from  3:30  to  7:15  hrs.  UT 
on  day  165.  The  radiometer  data  were  used  to  provide  the 
absolute  calibration  far  the  interferometer ,  but  as  the 
figure  shows  the  trends  in  the  interferometer  intensity 
curve  track  the  radiometer  curve  closely.  The  intensity 
peak  appears  at  about  5:dO  hrs.  UT  on  this  day.  The 
standard  deviation  of  the  interferometer  intensity 
calculation  is  about  2%  (the  standard  deviation  curve  is 
included  in  Appendix  C).  The  interferometer  was  realigned 
at  6:30  hrs.  UT .  Just  prior  to  real ignment,  the  standard 
deviation  of  the  intensity  increased.  After  alignment  the 
intensity  appeared  to  have  increased.  These  observations 
make  the  rapid  intensity  decline  at  6 : OO  hrs.  UT  suspect,  as 
an  instrument  alignment  drift  problem. 

Figure  5-9  presents  the  interferometer  OH  <3,1)  band 
intensity  for  UT  day  166.  Again  It  compares  favorably  with 
the  radiometer  curve  for  the  same  day  in  Figure  5-7.  The 
intensity  peak  occurs  on  this  day  at  5:30  hrs.  UT .  The  same 
reservations  about  alignment  drift  occurred  on  this  day  at 
6:00  UT  as  happened  on  the  previous  day. 
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The  calculated  i nterf erometer  OH  radiance  and 
temperature  curves  tor  the  other  OH  Meinel  bands  measured  in 
this  study  are  included  in  Appendix  C.  The  absolute  scales 
o-t  radiance  -for  the  other  OH  bands  were  based  on  the  (3,  1) 
band  (the  radiometer  was  calibrated  in  terms  of  the  (3,1) 
band)  and  using  the  band  intensity  ratios  developed  by 
Turnbull  and  Lowe  [19831  ( I (4, 2) SI (3,  1 ) x 1 . 03 , 
I  (8,5) =1  (3,  1  > xO.  14,  I  (7,41=1  (3,  11x0.09)  .  Only  when  the 
i  nter  f  erometer  is  viewing  in  the  zenith  are  the  intensity 
curves  plotted  in  absolute  terms.  The  geometry  o-f  looking 
at  low  elevation  angles  -tor  the  other  data  precluded 
inclusion  o-f  absolute  scales.  Also  shown  in  the  Appendix 
are  the  standard  deviations  -for  all  the  calculations. 

Rotational  Temperatures  for 
Interf erometer  Zenith  Data 

Included  in  Figures  5-8  and  5-9  are  the  calculated 
rotational  temperatures  for  days  165  and  166  when  the 
interferometer  is  viewing  in  the  zenith.  The  temperature 
tor  day  165  indicates  a  rise  in  temperature  and  intensity, 
beginning  at  about  4:00  hrs.  UT ,  with  a  mean  temperature  of 
about  175  *K.  The  total  rise  in  temperature  is  15  *K  with 
the  peak  occurring  at  5:20  hrs.  UT.  The  temperature  curve 
for  day  166  shows  the  same  increasing  trend  as  on  the 
previous  day.  The  mean  is  about  170  »K  with  a  total  rise  in 
temperature  of  15  *K.  The  temperature  peaks  on  day  166  at 
5:10  hrs.  UT .  On  both  days  as  the  interferometer  views  the 
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zenith,  the  peak  of  the  temperature  curve  precedes  the  peak 
of  the  intensity  curve  by  about  20  minutes  (discussed  more 
1 n  Chapter  VI ) . 


Rotational  Temperature 
Smoothing  Algorithm 

The  curve  of  calculated  rotational  temperatures  appeared 
quite  noisy  at  times.  Consequently,  smoothing  algorithms 
were  employed.  The  smoothing  technique  used  is  a  3  frame 
wide  sliding  window  where  the  averaging  is  a  weighted  one. 
Only  3  data  points  were  used  so  as  to  not  degrade  the 
temporal  resolution  of  the  data  (this  is  about  a  2  minute 
time  window)  and  still  provide  sufficient  smoothing.  The 
weighting  is  accomplished  using  the  reciprocal  of  the 
standard  deviation  squared  (or  reciprocal  of  variance)  for 
each  data  paint  as  its  weighting  factor  in  the  running  sum! 

therefore,  if  a  particular  data  point  has  a  large 

uncertainty  it  is  weighted  less  in  the  average.  This 

algorithm  was  chosen  because  in  data  sets  with  large 
differences  from  data  point  to  data  point,  with  significant 
differences  in  the  variance  of  each  point,  this  weighting 
technique  best  identifies  the  mean  curve  through  all  data 
points  IBevington  19693. 

The  specific  algorithm  used  for  the  rotational 

temperature  smoothing  is 
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where 

T  =  weighted  average  temperature  at  time  i, 

T  =  temperature  at  time  i, 
i 

c  =  standard  deviation  of  temperature  T  . 

T  i  1 

The  standard  deviation  on  the  smoothed  curve  is  also 
changed  due  to  the  weighting  and  is  recalculated  using 
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where 


=  averaged  standard  deviation  at  time  i. 

Figure  5-10  shows  temperature  and  standard  deviation 
curves  before  smoothing  and  Figure  5-11  shows  the  same 
curves  after  smoothing.  All  the  temperature  curves  in 
Appendix  C  show  a  set  of  temperature  curves  for  both  before 
and  after  smoothing. 

Interferometer  Recorded  Structure 

This  research  project  was  undertaken  in  an  attempt  to 
quantify  the  Intensity  modulations  and  anticipated  OH 
rotational  temperature  modulations  associated  with  airglow 
structure  events  typified  in  Figure  1-4.  The  presentation 
in  this  section  is  the  result  of  the  processed 
i nter f erometer  OH  (3,1)  band  records,  which  were  observed 
during  the  airglow  structure  event  of  June  15,  1963.  The 
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Figure  3-10.  OH  <4, 2)  band  rotational  temperature 
standard  deviation  before  smoothing. 
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OH  (3,1)  band  data  have  the  highest  signal -to-noi se  ratio 
and,  therefore,  are  felt  to  be  the  most  reliable.  The 
presentation  is  given  in  three  segments  cor  respond inq  to 
when  the  instrument  was  moved  in  viewing  elevation  or 
azimuth.  A  complete  catalog  of  all  the  observed  OH  band 
intensities  and  rotational  temperatures  is  contained  in 
Appendix  C. 

The  inter-ferometer-isocon  camera  system  was  lowered  to  a 
viewing  elevation  angle  o-f  17*  and  an  azimuth  angle  of  328* 
at  about  7:30  hrs.  UT  (after  moonset)  on  this  date.  Figure 
show  the  intensity  modulations  (and  standard  deviation) 
as  seen  by  the  i nterf erometer  until  about  8:30  hrs.  UT .  The 
bright  band,  dark  band,  bright  band  sequence  show  a 
modulation  in  intensity  of  about  20%.  This  same  sequence  is 
also  shown  in  the  isocon  video  frames  in  Figures  5-2  through 
5-4.  Following  the  last  bright  band,  the  intensity  falls 
off  by  40%,  indicating  a  relatively  dark  band.  The 

simultaneous  plot  of  the  OH  rotational  temperature  shown  in 
Figure  5-13  has  a  mean  temperature  of  about  165  »K.  The 

modulations  seen  are  in  phase  and  correspond  to  the 

increases  and  decreases  in  intensity  in  Figure  5-12.  The 

magnitude  of  the  temperature  modulations  are  from  5  to  8  *K. 
The  standard  deviation  on  the  temperature  calculation  is 
about  ±3  • K . 

The  i nter f erometer  was  moved  to  a  different  viewing 
location  at  about  8:30  hrs.  UT  and  remained  there  until 

about  9:15  hrs.  UT.  The  new  viewing  angles  were  15.5* 
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Figure  5-12.  OH  (7,1)  band  relative  intensity  and  standard 
deviation,  viewing  angle  *  17*  El.  328*  A  z .  ,  day  1 66 ,  7 : 30 - 
6:30  hr  s .  UT. 
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Figure  3-13.  OH  (3,1)  band  smoothed  rotation  temperature 
and  standard  deviation,  viewing  angle  “  17*  El.  328*  A2.  , 
day  166,  7:30-8:30  hrs.  UT . 
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elevation  and  340*  azimuth.  The  intensity  modulations  at 
this  viewing  position  are  shown  in  Figure  5-14.  The  -first 
dark  band,  bright  band,  dark  band  sequence  in  this  -figure 
show  a  modulation  o-f  40%.  The  intensity  modulations  shown 
after  the  first  bright  band  are  much  smaller,  being  on  the 
order  of  10%.  The  video  frames  in  Figures  5-5  and  5-6 
correspond  to  the  first  dark  band  and  occur  just  after  the 
peak  of  the  very  bright  band  of  Figure  5-14.  The  rotational 
temperature  plot  for  this  same  time  period  is  shown  in 
Figure  5-15.  The  mean  temperature  is  about  165  *K.  with  a 
maximum  modulation  of  10  *K.  The  uncertainty  on  the 
temperature  calculations  is  about  ±4  *K  during  this  time 
frame.  The  temperature  and  intensity  modulations,  at  the 
beginning  of  this  measurement  period,  are  again  in  phase 
with  each  other.  The  10%  intensity  modulations  shown  in 
Figure  5-14  occurring  after  8:50  hrs.  UT,  however,  do  not 
show  any  discernible  modulations  in  the  temperature. 

The  last  observations  of  the  night  began  at  about  9:15 
and  lasted  until  10:15  hrs.  UT.  The  viewing  position  of  the 
interferometer  was  again  changed  to  15.5*  elevation  and  309* 
azimuth.  Figure  5-16  presents  the  intensity  modulation 
record  for  this  period.  The  curve  shows  the  first  sequence 
of  structure  exhibiting  a  modulation  of  40%  and  three  more 
bright  and  associated  dark  bands  with  modulations  of  about 
20%.  The  corresponding  rotational  temperature  plot,  Figure 
5-17,  once  again  has  a  mean  temperature  of  165  *K.  The 


temperature  modulations  are  from  5  to  8  *K  with 


each  rise 


STD  DEV  INTENSITY  C»0  RELATIVE  INTENSITY 
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Figure  5-14.  OH  <3,1)  band  relative  intensity  and  standard 
deviation,  viewing  angle  =  13.5*  El.  340»Az.,  day  166,  8:30- 
9: IS  hrs .  UT. 
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Figure  5-13.  OH  (3,1)  band  smoothed  rotation  temperature 
and  standard  deviation,  viewing  angle  =  15.3*  El.  340*  Az .  , 
day  166,  8:30-9:13  hrs.  UT. 


105 


In  temperature  having  a  cor respond i ng  rise  in  intensity  of 
at  least  20%.  The  uncertainty  on  the  temperature 
calculations  is  again  about  id  • K . 

The  presentation  in  this  chapter  of  the  interierometer- 
recorded  structure  tor  June  15,  1983  is  focused  on  the 
information  extracted  from  the  OH  <3,1)  band.  Examination 
of  the  Appendix  C  records  tor  the  other  OH  bands,  presents 
an  additional  observation.  The  two  4v-2  band  temperatures 
track  each  other  within  the  standard  deviation  of  the 
calculations.  The  two  av=3  band  temperatures  track  each 
other  within  the  standard  deviation  of  the  calculations. 
The  av=3  band  temperatures ,  however,  are  consistently  from 
15  to  32  • K  hatter  than  the  4v= 2  band  temperatures 
(discussed  more  in  Chapter  V/I  >  .  Figure  5-18  is  an  example 
with  additional  data  available  in  Appendix  C.  Table  5-1 
outlines  a  summary  of  the  results  presented  in  this  chapter. 

TABLE  5-1.  Summary  of  OH  airglow  structure  measurement 
results  far  June  15,  1983. 


1.  Apparent  structure  period  .  1  d±  1  minutes 

2.  Apparent  structure  temporal  wavelength  .  2d± 1  km 


3.  Apparent  structure  phase  velocity  .  .  .  28i 2  meters/sec 


d.  Intensity  modulations  measured . .  10  to  dO  % 

5.  Rotational  temperature 

modulations  measured . 5  to  10  *K 


6.  Phase  relationship  of  recorded  rotational 
temperature  and  intensity  measurements  ....  In  phase 

7.  Mean  zenith  OH  <3,  1)  band  intensity . 75  kR 
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Figure  5-17.  OH  (3,1)  band  smoothed  rotation  temperature 
and  standard  deviation,  viewing  angle  =  15.5*  El.  309*.  Az .  , 
day  166,  9:15-10:15  hrs.  UT . 
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Figure  5-18.  OH  (7,4)  band  smoothed  rotation  temperature 
and  standard  deviation,  viewing  angle  ■  15.5*  El.  340*  Az . , 
day  166,  8:30-9:13  hrs.  UT. 
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CHAPTER  VI 

DISCUSSION  OF  RESULTS 

Rotational  Temperatures 

The  mean  mesopause  temperature  -for  the  month  o-f  June  at 
a  mid-latitude  site  is  expected  to  be  about  170  *K  with 
variations  o-f  ±20  *K  possible  during  the  month  [NOAA  1976 
and  references  therein!.  Noxon  [1978]  also  recorded  OH 
Meinel  rotational  temperatures  at  Fritz  Peak,  Colorado  (idO* 
N)  during  May  1977..  During  the  last  tew  days  of  May,  he 
recorded  nightly  mean  temperatures  of  about  160  *K. 

The  mean  OH  rotation  temperatures  presented  in  Chapter  V 
are  for  an  observing  site  at  *32*  N  and  are  between  165  »K 
and  175  *K.  The  references  cited  above  suggest  that  these 
rotational  temperatures  are  typical  of  the  mesopause 
temperatures  expected  at  mid-latitudes  during  the  summer 
season . 

Examination  of  the  standard  deviation  plots  on 
temperature  (Chapter  V  and  Appendix  C)  reveals  typical 
values  in  the  range  3-7  • K.  This  uncertainty  is  nearly  as 
large  as  many  of  the  temperature  changes  obtained  from  the 
structure  measurements .  It  is  felt,  however,  that  much  of 
the  computed  standard  deviation  may  be  systematic  rather 
than  statistical.  The  model  used  for  the  determination  of 
rotational  temperature  is  based  on  the  assumption  that  OH 


m 


no 


rotational  populations  are  in  true  thermal  equilibrium  and 
are  thus  strictly  Boltzmann  distributed.  A  slight  deviation 
■from  this  assumption  would  cause  a  systematic  error. 
Another  possible  source  of  error  results  -from  the  assumption 
that  all  OH  airglow  radiation  is  emitted  -from  a  thin  uni-form 
layer;  whereas  in  reality,  the  layer  is  about  7  km  in 
thickness.  In  addition,  at  low  viewing  elevation  angles 
%10*  the  layer  geometry  is  much  more  complicated.  At  these 
low  elevation  angles  atmospheric  extinction,  van  Rhijn 
effect ,  and  curved  spherical  geometry  potentially  have  a 
significant  impact  on  the  interpretation  o-f  the 
measurements. 

The  interferometer  spectral  response  calibration  is  very 
sensitive  to  alignment.  The  instrument  typically  remained 
in  alignment  Tor  about  2  hours.  As  can  be  seen  -from  the 
increase  in  the  standard  deviation  as  a  function  of  time 
(7:30  to  10: 15  hrs.  UT,  day  166),  the  alignment  changed 
significantly  and  this  change  could  account  for  a  portion  of 
the  uncertainty.  Therefore,  the  temperature  modulations 
obtained  from  the  spectral  data  are  felt  to  be  more 
accurately  defined  than  is  suggested  by  the  standard 
dev i at i on . 

The  assumption  that  the  low  observation  elevation  angles 
associated  with  the  structure  measurements  did  not  unduly 
impact  the  computed  temperature  is  supported  by  the  fact 
that  the  temperature  at  6:45  hrs.  UT  on  day  166  was  163  •  K 
(see  Figure  5-9)  and  the  temperature  at  7:30  hrs.  on  the 
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same  day  was  166  *K  (see  Figure  5-13) .  These  temperatures 
represent  the  values  computed  just  before  and  just  after  the 
Interferometer  "look  direction*  was  changed  from  the  zenith 
to  near  the  horizon. 

Rotational  Temperature  and 
Intensity  Modulations 

The  ranges  of  temperature  and  intensity  modulations 
observed  in  the  OH  Meinel  airglow  structures  are  given  in 
Chapter  V.  The  "adiabatic  oscillation"  and  "IGW  modeling 
of  the  OH  Meinel  airglow  variations  mentioned  in  Chapter  1 
utilize  a  parameter  which  is  readily  calculated  -from  the 
intensity  and  temperature  modulations.  This  parameter  is 
the  ratio  o-f  the  change  in  emission  intensity  normalized  b/ 
the  mean  emission  intensity,  divided  by  the  change  in 
temperature  normalized  by  the  mean  temperature  and  is  deemed 
useful  in  studies  of  the  OH  airglow  structure  phenomena. 
The  parameter  is  usually  represented  by  the  Greek  letter  eta 
< «1  )  and  is  defined  as  follows: 

H  =  [(6l/I)/(4T/T) 1  ,  (6.1) 

where 

Al  *  change  or  modulation  in  emission  intensity, 

I  ”  mean  value  of  the  emission  intensity, 

6T  “  change  or  modulation  in  rotational  temperature, 

T  “  mean  value  of  the  rotational  temperature. 


The  value  of  »l 


Is  potentially  useful 


in  distinguishing 


between  chemical  processes  which  give  rise  to  the  OH  airglow 

emission  and  temperature  structure.  The  physics  o-f  this 

parameter  is  discussed  by  Krassovsky  119721  and  Ueinstock 

119781.  Pendleton  1  19851  has  summarized  the  essential 

■features  of  this  parameter  in  Figure  6-1.  In  this  figure, 

the  *1  value  is  plotted  versus  the  ratio  (H/H  ),  where  H  is 

the  appropriate  atmospheric  scale  height,  and  H  is  the 

scale  height  (near  85  Urn)  of  minor  species  "  x  “  .  Here  the 

letter  "  x  "  represents  either  oxygen  (0)  or  hydrogen  (H). 

The  simple  ad i abat i c -osc i  1  1  at i on  model  of  Krassovsky  119721 

yields  *1  values  which  are  independent  of  (H/H  )  ,  whereas  the 

gravity-wave  model  of  Ueinstock  119781  yields  (H/H  ) 

dependent  values.  The  range  of  tl  values  expected  on  the 

basis  of  values  of  (H/H  )  inferred  from  several  measured 

x 

atomic  oxygen  profiles  is  also  shown  in  the  figure.  The 
information  in  Figure  6-1  indicates  that  values  of  »1  in  the 
range  from  3  to  6  might  be  expected  based  on  current 
gravity-wave  modeling  and  the  ozone  hydration  process. 

Using  the  numbers  for  the  intensity  and  temperature 
modulations  and  means,  presented  in  Chapter  V,  the  range  of 
calculated  values  for  n  are  from  8  to  12.  These  values  are 
about  a  factor  of  2  greater  than  those  shown  in  Figure  6-1. 
In  view  of  the  relatively  large  standard  deviations  on  the 
temperature  determinations,  the  nominal  factor  of  two 
disparity  between  predicted  and  measured  f)  values  is  not 
deemed  significant. 
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Figure  6-1.  Expected  values  <or  q  (vertical  axis)  based  on 
gravity-wave  models  and  oxygen  measurements  [Pendleton 
19851  . 
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The  »l  values  calculated  tor  the  time  period*  when  the 
interferometer  was  viewing  In  the  zenith  were  obtained  bv 
assuming  the  pre-midnight  increase  in  both  temperature  and 
intensity  reflected  a  wave-liUe  disturbance.  The  zenith¬ 
viewing  *1  values  associated  with  the  major  pre-midnight 
<I,T)  fluctuations  on  UT  days  163  and  166  were  found  to  be 
consistent  with  «l  values  deduced  from  the  1  ow-e  1  evat  i  on  - 
angle  data.  The  consistency  of  these  two  sets  of 
calculations  lends  credence  to  the  idea,  once  again,  that 
viewing  near  the  horizon  had  little  impact  on  the  *1 
determinations  although  the  modulation  in  both  intensity  and 
temperature  may  have  been  impacted  by  the  geometry  of  the 
measurements . 

The  field  of  view  1FOUI  of  the  interferometer  is  0.8* 
full  field.  Consideration  of  the  OH  structures  in  the 
nominally  13*  by  15*  video  frames  suggests  that  the  0.8*  FOV 
of  the  interferometer  results  in  a  horizontal  spatial 
integration  over  about  one  half  cycle  in  the  quasi-period 
structures.  This  integration  will  degrade  the  horizontal 
spatial  resolution  of  the  interferometer  measurements.  If 
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about  1.3.  The  horizontal  Intensity  structure  is  expected 
to  be  more  complex  than  the  temperature  structure  tWeinstocU 
19781.  However,  if  a  similar  ‘degradation  factor*  were 


applicable  to  the  intensity  measurements,  then  the 
satisfactory  agreement  of  the  1  values  deduced  from  the 
zenith-viexinq  and  near - hor l z on  v 1 ew 1 ng  measurements  could 
be  exp  1  a  1 ned . 

Temperature  Differences  Observed 
Between  Bands 

Rotational  temperature  differences  between  the  hiqh-v' 
and  1 o« - v  *  levels  characterized  both  the  zenith-viewing  and 
near - hor i zon - v i ew i ng  measur emen t s .  The  high-v'  rotational 
temperatures  were  consistently  higher  than  those  obtained 
from  low-v'  bands.  The  bands  used  tor  detailed  comparison 
are  the  high  v'  (7,4)  band  and  the  low  v*  (3,1)  band.  These 
bands  were  chosen  because  of  favorable  instrument  response 
and  alignment  stability  In  the  spectral  regions  of 
occurrence.  The  t emper a t u r e  differences  observed  ranged 
between  IS  and  32  *K.  The  smaller  difference  applied  when 
looking  In  the  zenith,  and  the  difference  gradually 
Increased  throughout  the  night  as  the  telescope  was 
physically  moved,  affecting  the  alignment.  The  increase  in 
the  temperature  difference  is  largely  attributed  to 
instrument  alignment  drift  because  the  standard  deviations 
on  the  temperature  calculations  (which  also  reflect 
m 1 sa 1 i gnmen t )  increase  simultaneously  with  (and  at  about  the 
same  rate)  the  Increase  in  temperature  difference. 

Explanations  were  sought  for  this  difference.  It  was 
discovered  that  art  error  In  the  1  Inc  strength  constant  for 


the  OH  (7,4),  Pi (4)  line  had  been  entered  into  the 
processing  system.  The  constant  was  about  20*4  low  in  value. 
In  order  to  assess  the  impact  o-f  this  error,  a  synthetic  OH 
spectrum  was  generated  assuming  a  typical  temperature  with 
the  associated  error  entered,  and  a  Boltzmann  plot  was  made. 
A  1  ine  was  -fit  to  the  points  in  a  1  eas  t  -  squares  sense  and  an 
associated  temperature  extracted.  The  OH  (7,4)  band  model 
used  12  lines  in  the  -fit;  there-f ore,  because  the  model  fit 
is  1  east -squares  in  nature  the  error  from  this  incorrect 
constant  was  found  to  be  less  than  1%. 

Another  possible  explanation  of  the  high-v'  ,  low-v’ 
t  emperature  difference  is  slightly  different  (3:  1-3  Urn) 
emitting  altitudes.  The  ratio  of  AT  to  the  mean  temperature 
T,  if  different  for  the  high-v'  low-v'  measurements,  would 
support  such  a  difference  in  mean  emitting  heights 
[Pendleton  19851.  The  4T/T  ratio  within  each  observational 
time  frame  was  calculated  and  the  difference  in  the  ratio, 
between  the  high  and  low  rotational  levels,  was  found  to  be 
less  than  10%,  with  less  than  30%  difference  among  all 
frames.  These  differences  in  the  fiT/T  ratio  between  high-v' 
and  low-v'  bands  when  viewed  in  terms  of  the  calculated 
standard  deviations  does  not  provide  evidence  for 
differences  in  mean  emitting  altitude. 

Within  the  standard  deviation  of  the  calculations,  the 
rotational  temperature  differences  between  the  high-v'  and 
low-v'  levels  appears  to  be  real.  The  magnitude  for  this 
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d  i  f  f  erence ,  as  was  mentioned  in  Chapter  1,  is  within  the 
range  reported  by  Krassovsky  and  Shagaev  119771. 

Temperature  and  Intensity 
Phase  Relationship 

The  modeling  of  IGW’s  mentioned  earlier  predicts  that 
the  changes  in  the  OH  Meinel  rotational  temperature  should 
be  in  phase  with  the  IGW  [Hines  19603.  The  change  in 
Intensity,  however,  should  be  related  to  the  IGW  temporal 
structure  in  a  potentially  more  complicated  manner.  This 
situation  arises  partly  because  o-f  the  finite  chemical  time 
constant  associated  with  mesospheric  ozone.  The  appropriate 
time  constant  -for  the  cool  (T  a  160  »K)  summer -mesopause 
conditions  is  about  25  minutes  IPendleton  19853.  For  IGW 
periods  somewhat  in  excess  o-f  this  value,  chemical 
conversion  of  O  is  expected  to  be  a  significant  factor  in 

w 

the  phenomenology,  whereas  for  much  smaller  IGW  periods 
chemical  conversion  should  be  unimportant.  Examination  of 
Figures  5-8  and  5-9  shows  that  the  temperature  maxima  lead 
the  intensity  maxima  by  about  20  minutes.  This  apparent 
phase  difference  may  relate  to  the  aforementioned  time 
constant,  but  it  would  be  premature  to  draw  this  conclusion. 
It  is  suggested  that  additional  attention  be  given  to  this 
Interesting  possibility. 

The  small-scale  structures  observed  near  the  horizon  on 
UT  day  166  exhibit  an  in-phase  relationship  between 
intensity  and  temperature.  The  only  exception  to  this  is 


when  the  measured  intensity  modulation  was  less  than  about 
20%.  Under  these  conditions,  no  direct  correlation  between 
intensity  and  temperature  could  be  drawn.  Since  the 
measured  period  of  the  small-scale  OH  Meinel  structures  was 
signi-f  leant  ly  less  than  the  nominal  85-km  0^  time  constant, 
it  does  not  appear  that  the  in-phase  behavior  of  the  larger- 
amplitude  -fluctuations  is  necessarily  inconsistent  with  the 
zenith  measurements. 
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CHAPTER  Mil 

CONCLUSIONS  AND  RECOMMENDATIONS 

Over  v 1 ew 

The  goal  of  this  study  was  to  design,  develop,  and 
operate  an  instrument  system  capable  of  performing 
simultaneous  spatial,  spectral,  and  temporal  high-resolution 
OH  airglaH  measurements.  The  design  herein,  and  the 
resulting  data  demonstrate  the  ef  f  ec  t  i  veness  of  the 
technique.  An  airglaw  structure  event  which  occurred  on 
June  13,  1983  was  measured  with  the  interferometer- 

spectrometer  system.  OH  Melnel  intensities  and  rotational 
temperatures  were  obtained  for  the  peaks  and  troughs  of  this 
wavelike  structure. 

Cone  1  us  1 ons 

The  following  are  the  specific  accomp 1 i shments  of  this 
study.  The  areas  addressed  pertain  both  to  the  instrument 
designed  for  airglow  structure  measurements  and  to  the  data 
processing  techniques  used. 

1.  An  opt i ca 1 1 y-compensated  i nterferometer  for  high 

2 

throughput  <Aft*0.285  cm  sr ) ,  was  matched  to  a  large 
area  collector  iSO-cm  diameter)  to  narrow  the  field  of 
view  (0.8*).  A  noise  equivalent  spectral  radiance 
(NESR)  (sensitivity)  of  16  R/cm  1  at  1.3  sm  was 
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achieved.  In  comparison,  a  conventional  Michelson 
1 n ter  Aerometer -spec trometer  with  the  same  detector, 
collector  area,  resolution,  and  scan  time  would  have  a 
NESR  of  128  R/cm  a  factor  of  8  less  sensitive  than 
the  one  developed  for  this  study.  When  compared  with  a 
conventional  Ebert  spectrometer,  using  the  same  detector 
and  operated  at  the  same  resolution,  a  sensitivity  of 
208  R/cm  1  could  be  achieved,  a  factor  of  13  less 
sens i t i ve . 

2.  A  spectral  resolution  of  2  cm  1  was  sufficient  to 
resolve  the  OH  emission  line  structure  for  the 
extraction  of  OH  rotational  temperatures.  Based  upon  a 
rotational  line  separation  of  10  cm  1  and  the  Hammir  t 
apodization  function  used,  the  spectral  resolution  of 
the  Instrument  could  be  lowered  to  no  more  than  4  cm  1. 

3.  Based  upon  the  video  records,  the  bright  or  dark  bands 

of  OH  structure  subtend  about  1*  of  arc  at  these  low 
elevation  angles  <*  15* > .  The  apparent  temporal 

wavelength  was  24± 1  km,  with  a  period  of  14±1  minutes, 
and  an  apparent  phase  velocity  of  2812  meter s /second . 

4.  The  i nter f erometer  system  field  of  view  was  measured  at 
0.8*.  The  i nter f erometer  FOV  is  sufficiently  narrow  to 
independently  view  a  "bright"  or  a  "dark"  structure 
band.  Based  on  this  limited  data  set  and  simple  wave 
geometry,  the  field  of  view  could  be  as  large  as  7*  if 


these  structures  were  viewed  in  the  zenith. 
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5.  The  measured  intensity  modulations  (contrast  ratios)  -for 
the  OH  airglow  structures  ranged  between  20  and  40%  with 
i n ter f erometer  recorded  periods  of  14  minutes  The 
calculated  standard  deviation  was  typically  3%. 

6.  The  mean  calculated  OH  Meinel  rotational  temperature  -for 

the  af or emen t i aned  event  was  165  *K.  The  measured 

modulations  in  rotational  temperatures  associated  with 
the  changes  in  intensity  ranged  -from  5-10  *K  and  are  in 
phase  with  the  intensity  modulations.  Typical  standard 
deviations  on  the  rotational  temperature  calculations 
ranged  -from  2-7  *K.  The  mean  temperature  and  magnitude 
o-f  the  temperature  -fluctuations  are  consistent  with  both 
I GW  theory  and  previous  mesospheric  temperature 
measur emen  ts . 

7.  The  least-squares  model  used  to  extract  band  intensity 
and  rotational  temperature  provides  a  computational 
efficient  way  (convergence  to  final  values  occurred 
within  4  iterations  of  the  fitting  routine)  to 
simultaneously  derive  these  values.  The  model  also 
provides  a  measure  (standard  deviation)  of  how  well  the 
data  fit  a  Boltzmann  distribution. 

The  insight  needed  to  identify  what  is  being  observed 
from  the  airglow  layer  can  only  be  provided  by  the 
camera  (or  similar)  video  system.  The  measurement  of  OH 
airglow  structure  events  with  the  interferometer  system 
would  not  be  possible  without  the  simultaneous  use  of 
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the  lsocon  camera,  because  with  the  interferometer 
alone,  exactly  what  was  being  viewed  would  be  unknown. 

Recommendations  -for  Future  Research 

The  interferometer  worked  as  designed  and  proved  to  be 
an  excellent  tool  -for  this  type  of  study.  The  model 
developed  for  the  extraction  of  rotational  temperatures  and 
intensities  is  an  accurate  technique  and  provides 
computational  flexibility.  However,  several  suggestions  are 
made  for  consideration  for  future  work. 

1.  The  optical  path  within  the  i nterf erometer  is  very 

complex.  There  are  20  optical  surfaces  through  which 
the  incoming  energy  must  pass  before  reaching  the 
detector.  Assuming  a  typical  loss  of  AX  per  surface, 
80%  of  the  incoming  signal  is  lost  before  reaching  the 
detector.  This  complex  optical  path  should  be 

redesigned  to  minimize  the  number  of  optical  elements 
and  optically  coat  the  remaining  elements  to  minimize 
reflection  loss. 

2.  The  physical  size  of  the  instrument  should  be  reduced  to 
facilitate  portability  to  remote  sites. 

3.  The  instrument  is  very  sensitive  to  optical  alignment. 
The  optical  components  need  to  be  mounted  in  a  more 
stable  manner  to  hold  their  location  better. 

4.  The  alignment  is  very  sensitive  to  temperature.  It  is 
necessary  therefore  to  temperature  control  the 
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interferometer  environment  to  preclude  instrument 

changes  as  the  outside  temperature  fluctuates. 

5.  The  telescope  pointing  system  needs  to  be  automated. 
The  sensitivity  of  alignment  to  physical  motion  as  well 
as  the  need  to  point  the  instrument  to  an  area  of 
interest  demand  that  the  positioning  of  the  system  be 
automated . 

6.  An  infrared  camera  system,  like  the  one  provided  by  the 

University  of  Southampton,  needs  to  be  permanently 

incorporated  as  part  of  the  i nterf erometer  system.  The 
interferometer  cannot  effectively  gather  data  on  airglow 
structure  if  the  location  of  the  structure  is  not  known 
and  the  camera  provides  this  input. 

7.  The  instrument  must  be  provided  a  better  means  of 

calibration.  The  blackbody  sources  used  in  this  study 
give  a  reasonable  indication  of  alignment  and  instrument 
response  but  as  the  system  alignment  drifts  the 

calibration  is  less  meaningful.  Perhaps  a  technique 
utilizing  OH  spectral  line  pair  ratios  which  are 

independent  of  rotational  temperature  but  sensitive  to 
alignment  could  be  used  as  a  dynamic  measure  of 
instrument  alignment. 

8.  A  recommendat i on  is  made  to  investigate  other  detectors, 
in  order  to  extend  the  ability  of  the  interferometer, 
with  a  wider  spectral  bandwidth  and  higher  sensitivity. 
The  RCA  detector  used  in  this  study  is  an  excellent 


detector  where  it  is  sensitive  but  is  somewhat  limited 
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in  spectral  bandwidth.  A  larger  detector  could  also  be 
used  to  increase  the  throughput,  although  the  throughput 
is  now  almost  limited  by  the  size  of  the  interferometer 
opt i cs . 

9.  A  new  data  system  should  be  developed  to  record  the 
interferometer  data  allowing  for  at  least  the 
digitization  o-f  the  data  during  recording.  The  analog 
tapes  are  bulky  and  expensive  but  more  importantly, 
playing  back  the  tapes  -for  data  reduction  is  too  time 
consum i ng . 

10.  The  calibration  curves  provided  to  the  model  derived 
from  occasional  alignment  processes  with  a  tungsten 
b 1 ackbody  could  be  improved.  As  was  mentioned  earlier, 
a  dynamic  calibration  using  in-formation  inherent  in  the 
spectrum  could  be  used  to  better  adjust  the  model  to  the 
instrument  response. 

11.  Examination  of  Figures  3-3  and  3-4  shows  that  the  phase 
correction  used  to  eliminate  the  chromatic  effects  of 
the  instrument  on  the  data  works  well  on  the  slowly 
varying  blackbody  curve  but  the  negative  information  on 
the  spectrum  makes  the  technique  suspect  when  applied  to 
the  rapidly  changing  data.  A  phase  correction  technique 
which  operates  in  the  i  nter f erogr am  "domain*  where  the 
shaping  could  be  done  with  simple  multiplication  could 
prove  to  be  more  accurate. 
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12.  The  extraction  of  line  amplitudes  from  the 

interferometer  could  be  improved  in  two  areas.  First, 
the  approximate  locations  of  the  spectral  lines  cf 
interest  are  found  in  a  manual  manner.  A  template  using 
a  synthetic  spectrum  could  be  designed  for  each  OH  band 
and  a  correlation  routine  could  be  used  to  automatically 
search  the  raw  data  far  the  location  of  the  lines. 
Secondly,  the  apodization  routine  used  to  extract  the 
actual  line  amplitude  from  the  data  should  be  modified 
to  calculate  a  line  area  rather  than  amplitude.  The 
area  routine  would  provide  for  the  averaging  out  of 
noise  whereas  the  amplitude  routine  always  searches  for 
the  most  positive  peak. 

13.  The  model  should  be  modified  to  to  use  the  "R"  branrhp? 

of  the  OH  bands.  To  do  this  the  molecular  constants  for 
both  and  would  need  to  be  averaged  as  one  because 

the  interferometer  does  not  resolve  the  two  groups. 
This  addition  should  add  more  accuracy  to  the  model 
because  the  Q  branches  are  the  largest  lines  within  each 
band  . 

Id.  The  model  could  be  modified  to  include  a  third  variable, 
water  column  content,  based  on  the  several  lines  within 
the  OH  bands  which  are  severely  affected  by  water 
absorption.  The  model  then  could  provide  additional 
information  about  the  atmosphere. 
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Append i *  A 
OH  T rans i t  i  ons 

This  appendix  briefly  describes  why  the  radiation  from 
the  OH  molecule  is  so  complex.  The  mar,/  spectral  lines 
generated  by  the  excited  radical  are  depicted  in  figure  *4  1. 

The  mo lecule  rotates  and  vibrates  simultaneously,  and  each 
of  the  mo  t  i  on  s  are  qu  antifled  arid  inter  ac  tut. 

The  total  molecular  angular  momentum  without  electronic 

■4 

spin  K  is  also  quantized  and  is  identified  by  the  quantum 

*4 

number  k.  The  quantity  h  is  comprised  ut  two  pci  ts 

■*  a  » 

K  -  N  ♦  A.  .  ui .  l 

-» 

The  vector  M  is  the  nuclear  angular  momentum  and  the  vector 

*4 

A  is  the  angular  momentum  of  the  orbiting  electron  cloud 
projected  onto  the  internuclear  axis.  The  quantum  number  A 
associated  with  the  electronic  orbital  momentum  can  take  on 
a  value  of  *1  or  -1  depending  upon  which  way  the  electron 
cloud  is  orbiting  with  respect  to  the  nuclear  rotation.  The 
double  degeneracy  of  A  leads  to  the  so-called  A  splitting  of 
each  state;  however,  the  split  is  less  than  1  cm  1  at  low 
rotational  speeds  [Baker  19781,  which  is  less  than  the 
instrument  resolution  used  for  this  study,  therefore  the  A 
split  lines  will  be  considered  as  one. 

The  quantum  number  K  can  take  on  values  k  -  1,2,3,... 

The  selection  rule;  however,  is 
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A  K  -  O  ,  ♦  1  .  ( A  .  2  ) 

f  r  nl  Ifctinri  cj  f  lir.es  within  each  band,  grouped  actor  ding 

to  their  respective  AK,  are  called  branches.  The  branch 

with  At  =  0  is  called  the  Q  branch,  that  -for  AK  =  +1  is  the 

P  branch,  and  that  (or  aK  =  -1  is  called  the  P  branch. 

The  OH  molecule  has  an  odd  number  of  electrons.  This 

unbalance  results  i  ri  a  net  electronic  spin  angular  momentum 
* 

P  and  is  represented  by  quantum  number  S.  The  odd  electron 
gives  rise  to  an  even  multiplicity  2S+1.  Since  the  total 
number  of  electrons  is  odd,  S  is  half  integral  (S  =±1/2), 
each  transition  state  is  a  doublet.  It  is  sometimes 
convenient  to  consider  the  total  electron  angular  momentum  A 
as  a  separate  entity.  The  total  electron  angular  momentum 
i  s 

A  =  A  «•  S  .  (A. 3) 

Therefore,  each  vibration-rotation  transition  will  split 
into  two  separate  spectral  lines  according  to  whether  A=3/2 
or  (1  =  1  t  2  . 

The  OH  molecule  is  very  light  and  as  a  consequence  the 
odd  electron  spin  3  is  only  weakly  coupled  to  the 
iriternuclnar  axis.  The  molecule  is  therefore,  modeled  as 
Hu  nd ’ s  case  (b)  CHertzberg  19711.  The  total  molecular 
ar.gu  1  ar  momentum  3  can  now  be  formed 


-» 

J 


a  -* 
K  +  S 


(  A.  4  ) 


V*1 

r 

■  V* 

v,'*i 
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As  can  be  seen  from  Equation  A.  4,  -tor  each  value  o-f  K  there 
are  two  values  -for  J.  Each  branch  o-f  the  OH  spectra  must 
also  take  on  two  values.  If  fl  = 3 / 2 ,  then 

J  =  K  +  1/2  =  1.5,  2.5,  3.5,  ...  .  (A. 5) 

These  values  for  J  lead  to  a  set  of  spectral  lines  known  as 

Pj,  Qj,  and  R  branches.  If  fl=l/2,  then 

J  =  K  -  1/2  =  0.5,  1.5,  2.5,  ...  .  (A. 6) 

These  values  for  J  give  rise  to  a  set  of  spectral  lines 

known  as  P  .  <3  ,  and  R  branches.  Additional  information  on 

2  2  2 

the  physics  of  OH  molecule  is  readily  available  in  the 
literature,  among  them  are  Baker  C1978J,  Hertzberg  119711, 
and  Mies  119743. 


APPENDIX  B. 


RCA  LIMITED  GERMANIUM  DETECTOR  SPECIFICATIONS 

RCA  supplied  the  following  specifications  for  their  limited  germanium 
detector: 

TE-71  Series 


NEP( 1 .42  ,10  ,1 )  =  10"14  WHz-1/2 
TM-71  Series 

NEP ( 1 . 42, 500,1)  =  10-14  WHz“1/2 
Area  =  0.20  cm2 
Cooled  to  77  K 

Responsivity  "  107  -  108  VW"1 


In  order  to  achieve  good  noise  performance,  it  is  necessary  to  cool  both  the 
detector  and  the  preamplifier. 

Special  Optics:  The  standard  window  is  quartz;  different  window  material, 
special  filter  or  condensing  optics  can  be  fitted  provided  no  substantial 
mechanical  redesign  is  necessary.  The  normal  field  of  view  for  the  detector 
is  close  to  90°. 

If  the  user  desires  to  operate  the  TM-71  system  at  very  high  frequencies,  it 
is  possible  to  trade  NEP  for  frequency  response.  Other  special  features  may 
be  added  to  these  developmental  systems. 
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The  responsivity  varies  according  to  the 
spectral  curve.  Some  variation  in  this 
curve  is  possible  according  to  the 
customer  requirements. 


Figure  B-l.  Typical  detector  spectral  response.  The  responsitivity  varies 
according  to  the  spectral  curve.  Some  variation  in  this  curve 
is  possible  according  to  the  customer  requirements. 


TABLE  B-l.  Germanium  Detector  Technical  Data 


TE  Series  71 


TM  Series  71 


Operating  Temperature 

Noise  Equivalent  Power 
(guaranteed),  WHz-1/2 

Best  NEP  (previously 
achieved) ,  WHz-  V2 

Noise  level  (at  output) 

Impedance  level  (at  output) 

Responsivity,  VW_1 

Linear  range  for  power,  W 

Useable  1 imit  (power) ,  W 

Frequency  caracteri Stic 

Detector  area  (circular),  cm2 

Typical  Operating  Volts,  V 


77  K 

NEP(1. 42,10, l)=lxl0-i- 
NEP(1.42,23,l)=lxl0"15 
~lyV 

-sooon 
~5xl0"7 
10" 14  -  10-10 
-6x10" 9 
f"l 

0.2 

-10 


77  K 

NEP(1.42,500,l)=lxl0~14 

NEP  ( 1.42,4000,1 )  =  1x10" 1 5 

~10PV 
~500ft 
~5xl0"7 
10" 1 4  -  10-10 
~2xl0-9 

flat,  3dB  at  500  Hz 

0.2 

-10,  +10 


1/4" 

VACUUM 

FITTING 


I 


Rl  < 


.>.V  v  V  -t 


Append i  x  C 


Interferomet«r  Data  Catalog 

The  interferometer  data  for  observations  taken  on  June 
13-15,  1983  are  contained  within  this  appendix  in  their 

entirety.  All  data  presented  were  recorded  at  Sacrament 
Peak,  New  Mexico.  The  figures  are  all  organised 

chronological ly  beginning  on  day  165  at  3:30  hrs.  UT  ani 
proceeding  through  day  166,  10:15  hrs.  UT . 

Both  of  the  observation  days'  records  begin  with  t  t> . - 
interferometer  viewing  in  the  zenith.  At  moonset  on  day  inn 
the  interferometer  and  camera  systems  were  lowered  to  view 
near  the  horizon.  The  low  elevation  viewing  period  dur  i n 
day  166  is  divided  into  three  time  frames,  corresponding  1 1. 
when  the  interferometer  was  adjusted  in  viewing  location. 
The  time  period  and  viewing  position  are  identified  in  earn 
figure  caption. 

Within  each  time  segment  of  the  data  presentation,  t he 
figures  are  organized  according  to  the  OH  Meinel  band  ♦ r  on 
which  the  data  were  calculated.  First  is  the  OH  ( *4 , 2 )  band, 
second  the  OH  <3,11,  third  the  OH  (8,5)  band,  and  fourth  the 
OH  (7,4)  band.  Each  of  the  band  groups  shows  curves  for 
first  the  intensity,  second  the  rotational  temperature,  and 
third  the  smoothed  rotational  temperature  calculation  wit), 
the  associated  standard  deviation  for  each. 
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Figure  C-15.  OH  (4,2)  band  smoothed  rotational  temperature 
and  standard  deviation,  viewing  angle  =  zenith,  day  166, 
3:30-6143  hrs .  UT. 


Figure  C-16.  OH  (3,1)  band  radiance  and  standard  deviation, 
viewing  angle  =  zenith,  day  166,  3:30-6:45  hrs.  UT. 
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Figure  C-17.  OH  (3,1)  band  rotational  temperature  and 
standard  deviation,  viewing  angle  =  zenith,  day  166,  3:30- 
6143  hrs.  UT. 
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Figure  C-21.  OH 
and  standard  di 
3:30-6:43  hrs.  U" 
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Figure  C-22.  OH  (7,4)  band  radiance  and  standard  deviation, 
viewing  angle  ■  zenith,  day  166,  3:30-6:45  hrs.  UT. 
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Figure  C-23.  OH  (7,4)  band  rotational  temperature  and 
standard  deviation,  viewing  angle  =  zenith,  day  166,  3:30- 
6 . 43  hrs.  UT. 
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Figure  C-24.  OH  (7,4)  band  smoothed  rotational  temperature 
and  standard  deviation,  viewing  angle  =  zenith,  day  166, 
3:30-6:43  hr s .  UT. 
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Figur#  C - 28 .  OH  (3,1)  band  rtl*t lvt  intensity  and  standard 
deviation,  viewing  angle  ■  17*  El.  328*  Az .  ,  day  166,  7:30- 
8:30  hr s .  UT. 
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F 1 gure  C-29.  OH  (3,1)  band 
standard  deviation,  viewing 
166,  7130-8:30  hr s .  UT. 
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Figure  C-30.  OH  (3,1)  band  smoothed  rotational  temperature 
and  standard  deviation,  viewing  angle  =  17*  El.  328*  Az.  , 
day  166,  7:30-8:30  hrs.  UT. 
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Figure  C-31.  OH  (8,3)  band  relative  intensity  and  standard 
deviation,  viewing  angle  *  17*  El.  328*  Az. ,  day  166,  7:30- 
8:30  hr s .  UT. 


230 


SACRAMENTO  PEAK 
HEU  MEXICO 


DAY  1661  JUNE  11-1S.  63 
OH  BAND  C  6- 5  ) 


TIHECUT) 


UJ 
0 c 
z> 


<E 

QC 

UJ 

a. 

c 

ui 


> 

UJ 

Q 

Q 

*- 

u> 


SACRAMENTO  PEAK  DAY  166«  JUNE  11-lS.  83 


TIME  CUT  3 


Figure  C-32.  OH  <8, 5)  band  rotational  temperature 
standard  deviation,  viewing  angle  =  17*  El.  328*  Az.  , 
166,  7:30-8:30  hrs.  UT . 
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Figure  C -34 .  OH  (7,4)  band  relative  Intensity  and  standard 
deviation,  viewing  angle  =*  17*  El.  328*  Az.,  day  166,  7:30- 
8:30  hrs.  UT . 
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Figure  C-35.  OH  <7, A)  band  rotational  temperature 
standard  deviation,  viewing  angle  “  17*  El.  328*  Az., 
166,  7:30-8:30  hr s .  UT. 
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Figure  C-37.  OH  (4,2)  band  relative  intenmity  and  standard 
deviation,  viewing  angle  *■  13.3*  El.  340*  Az.,  day  166, 
8:30-9: 13  hr*.  UT. 


v.i1}  I'ii'ii  c  ,^ms*vv — 

■  .;.V  .*  •  J  j'.  )  *  ’  i  *  ■  *  .  t 


;v  y 


V,  r.T, 


STD  DEV  TEMPERATURE  tK)  TEMPERATURE  IK) 


SACRAMENTO  PEAK  OAT  J66«  JUNE  11-1S.  83 


TlMEtUT) 


SACRAMENTO  PEAK  DAY  J66«  JUNE  11-1S.  83 


TIME(UT) 


Figure  C-38.  OH  (4,2)  band  rotational  temperature  and 
standard  deviation,  viewing  angle  =  13.5*  El.  340*  Az.  ,  day 
166,  8:30-9: 13  hrs.  UT . 
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Figure  C-40.  OH  (3,1)  bend  relative  intensity  end  standard 
deviation,  viewing  angle  »  19.9*  El.  340*  Az. ,  day  166, 
8:30-9: 19  hrs .  UT. 
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166,  8:30-9:13  hr • ,  UT . 
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Figure  C-44.  OH  (6,3)  band  rotational  temperature  and 
standard  deviation,  viewing  angle  “  13.3*  El.  340*  Az.,  day 
166,  6:30-9:13  hrs .  UT. 
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Figure  c-43.  OH  (8,3)  band  smoothed  rotational  temperature 
and  standard  deviation,  viewing  angle  “  13.3*  El.  340*  Ax., 
day  166,  8:30-9:13  hrs.  UT. 
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Figure  C-46.  OH  (7,4)  band  relative  intanaity  and  standard 
deviation,  viewing  angle  “  19. S*  El.  340*  Az . ,  day  166, 
6: 30-9: 13  hre.  UT. 
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Figure  c -47.  OH  <7, 4)  band  rotational  temperature  and 
standard  deviation,  viewing  angle  “  13.3*  El.  340*  Az.,  day 
166,  8:30-9: 13  Are.  UT. 
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Flgura  C-40 .  OH  <7, 4)  band  smoothad  rotational  tamparatura 
and  standard  davtation,  viawing  angla  ”  15.5*  El.  340*  Az., 
day  166,  8:30-9:15  hra.  UT. 
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Figure  C-49.  OH  (4,2)  bard  relative  intensity  and  standard 
deviation,  viewing  angle  ■  15.5*  El.  309*  Az.,  day  166, 
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Figure  C-30.  OH  (4,21  band  rotational  temperature 
standard  deviation,  viewing  angle  *  13.3*  El.  309*  Az 
166.  9:13-10:13  hrs.  UT. 
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Figure  C-51.  OH  (4,2)  band  smoothed  rotational  temperature 
and  standard  deviation,  viewing  angle  “  15. S*  El.  309*  Az., 
day  166,  9:13-10:13  hrs.  UT. 
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Figure  C-33.  OH  (3,1)  band  rotational 
standard  deviation,  viewing  angle  =  15.5* 
166,  9:15-10:13  hrs.  UT. 
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Figur*  C-54.  OH  (3t 1)  band  aaoothad  rotational  t»np«rttur* 
and  standard  deviation,  vlauing  angl*  -  19.9*  El.  309*  Ai.  , 
day  166,  9:19-10:19  hrs.  UT. 
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Pigur#  C-S6.  OH  <8,S)  band  rotational  tmptratur* 
standard  deviation,  viewing  angle  ■  13.3*  El.  309*  Az, 
1*A,  9:  13-10:  13  hr s .  t'T. 
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Figure  C-37.  OH  (8,5)  band  smoothed  rotational  temperature 
and  standard  deviation,  viewing  angle  »  15.5*  El.  309*  Az., 
day  166,  9:13-10:13  hrs.  UT. 
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Figure  C-38. 
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OH  17,4)  band  rtlativt  Intensity  and  standard 
viewing  angle  ”  IS. 5*  El.  309*  Az. ,  day  166, 
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Figure  c-39.  OH  (7,4)  band  rotational  temperatur< 
standard  deviation,  viewing  angle  *  13. S*  El. 

166,  9:13-10:13  hr s .  UT. 
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Figure  C-60.  OH  (7, A)  band  smoothed  rotational  temptratur* 
and  standard  deviation,  viewing  angla  -  19.9*  El.  309*  Az. , 
day  166,  9: 19-10:19  hr*.  UT. 


